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1. Introduction 


Background 


Automatic/Aided Target Recognition (ATR) systems are a developing class 
of machine vision devices that scan a field of regard (FOR), process that infor- 
mation, then report potential objects of interest (targets) to a human operator or 
to another automatic device. It is the intent of the U.S. Army to equip existing 
and next-generation attack and reconnaissance helicopters with ATR systems. 
Performance goals for ATR systems are that they have very high probability of 
locating and reporting valid targets and very low probability of reporting 
invalid or false targets for all conditions under which they will operate. 
Developing ATR systems use sensor information obtained from passive visible 
and thermal infrared imagers, millimeter wave (MMW) radar, and laser range 
finders, and non-sensor information such as digital map (i.e. terrain) data and 
location information from giobal positioning systems (GPS). 


Early testing of developing ATR systems has identified problems in achiev- 
ing required performance goals. Systems have exhibited inconsistent perfor- 
mance over their intended operational environment, low probability of 
detection, and high false alarm rates. Further, detection probability and false 
alarm rates have demonstrated a high degree of sensitivity to terrain and 
weather (collectively referred to as environmental) conditions, particularly for 
the passive sensors. Clearly, it is imperative to test under a broad range of 
environmental conditions during the development cycle. Testers and evalua- 
tors in the ATR development community are now burdened with determining 
the following: 


a. Which continental United States (CONUS) test sites should be used for 
testing and when they should be used. 


b. How to compare ATR performance results from different CONUS test 
sites—which sites represent more difficult conditions. 


c. Which CONUS test sites are most analogous to potential theater of 
operation sites. 
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d. How to specify environmental conditions in an ATR system perfor- 
mance test. 


The Environmental Characterization for Target Acquisition (ECTA) Pro- 
gram was initiated in an attempt to address these types of concerns, The pri- 
mary goals of the ECTA Program are to develop and apply a methodology to 
quantify "scene complexity” of available candidate ATR testing and training 
Sites and to establish methods for specifying environmental conditions for 
thermal and visible electro-optical (EO) imager systems and 35-GHz radar 
sensor systems. 


Approach 


Almost all ATR systems use pattern recognition techniques to detect targets 
within an imaged background area. This is true whether the system uses a 
passive EO imager or an active EO imager (e.g. laser) or radar sensors. The 
ATR logic filters the entire image for targetlike regions which it analyzes in 
greater detail to make first-leve] (detection) target acquisition decisions. 
Detection is the only stage in the target acquisition process that examines the 
entire image; all subsequent stages (classification, recognition, and identifica- 
tion) use only the regions of interest. Therefore, the ultimate success of target 
acquisition depends on the system’s ability to separate targets frorn background 
features that can have similar signature characteristics. 


Based on these considerations, an approach was formulated for processing 
imagery and signals in specific wavebands to determine the distribution of 
“targetlike" features within the background scene that could result in poor ATR 
performance (missed detections and false targets). Scenes having a high dens- 
ity of targetlike features are deemed to have high scene complexity and those 
with few such features are deemed to have low scene complexity. These 
image characterization measurements, or metrics, were determined using ther- 
mal and visible image data, and 35-GHz radar data representing a systematic 
sampling of commonly used CONUS test sites at different times of the day 
and of the year. These sites would be subjected to intensive measurement of 
key terrain and meteorological attributes during imaging periods. An analysis 
would be conducted to relate these measured attributes (referred to as ground- 
truth conditions) to values of the selected metrics for each specific wavelength 
and to devise a means of categorizing and predicting metric levels for these 
sites. Metrics are used here as indicators of the level of scene complexity, 
which is the relative density of targetlike features a background scene contains. 


Scope and Structure 


Several constraints were placed on the conduct of the program as follows: 
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a. Field measurements would be conducted at five sites commonly used 
for testing Army aviation ATR systems. 


b. Each site would be visited twice—once durin” leaf-on conditions and 
once during leaf-off conditions (for deciduous vegetated sites). 


c. Data collected during each site visit (excursion) would be limited to a 
3-day period that would include a single 24-hr imaging period. 


d. Existing U.S. Army Engineer Waterways Experiment Station (WES) 
instrumentation and imaging assets would be used (includes only visible 
and thermal imagers; no radar sensor was available), 


e. Imagery (35-GHz radar) collected by the Martin-Marietta Corp.'s 68D 
radar system during the Multi-Sensor Fusion Demonstrziion at Fort 
Hunter Liggett, CA, would be used, and radar analysis would be con- 
ducted as a stand-alone task, unrelated to the visible and thermal 
imaging tasks listed in items a and b above. 


The ECTA Program is documented in three separate reports. Report 1 
(Berry, Rivera, and Sabol 1993) describes ground-truth measurements made at 
each of the selected sites. That report describes procedures used for collection 
of terrain, meteorological, and radiometric data; appendixes contain a data 
summary for each site. The study reported herein (second in the series) 
describes analysis procedures for evaluating the relationships between the site’s 
ground-truth measurements and statistical characteristics of resulting imagery; 
results of these analyses are presented and discussed. The third report (Curtis 
and Sabol 1993) in this series describes the radar scene analyses conducted 
using the Martin-Marietta radar data collected at Fort Hunter Liggett, CA. 
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2 Methodology 


The overall methodology used in this project is illustrated in Figure 1. 
Field measurements and imagery were systematically collected from selected 
test sites (step 1). Detailed procedures for obtaining these data are described 
in ECTA Report 1 (Berry, Rivera, and Sabol 1993). A cursory description of 
these procedures is provided here for completeness. Image metrics indicative 
of scene complexity were selected and implemented in software (step 2). 
Imagery was processed to compute a set of metrics values for each individual 
image; these values were put into a database containing the corresponding 
weather and terrain attributes associated with each image (step 3). Finally, 
statistical analyses were performed to identify environmental and scenario vari- 
ables that most affect the scene complexity metrics values and to develop 
techniques for predicting scene complexity metrics levels from environmental 
and scenario data (step 4). Each of these steps is discussed in the following 
sections. 


Data Collection Procedures 


Data were collected at five CONUS test sites commonly used for Army 
aviation systems testing and evaluation. The sites were selected to allow data 
to be collected as an analog of temperate, coastal plains, and semiarid envi- 
ronmental conditions as indicated in Table 1. Temperate sites were sampled 
twice, once during the leaf-on period and once during the leaf-off period. 
Other sites were sampled only once since vegetation thereon does not change 
appreciably during the ycar. 


At each site, a single sensor location was selected to view a wide FOR and 
to achieve a viewing geometry approximating a helicopter scenario. This 
scenario is characterized by: (a) sensor position 5 to 20 m above the imaged 
terrain surface; (b) camera pointing angles ranging from horizontal to several 
degrees below horizontal; and (c) the distance of the center of the images 
ranging from several hundred meters to several kilometers. The selected FOR 
was divided into 11 to 22 contiguous 2.5°x2,5° "scenes" which were imaged by 
the suite of imagers over a 24-hr period at a 2-hr interval starting immediately 
before dawn on the second day of each excursion. 
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1. Systematically collect 
imagery and ground-truth 
data at selected sites 


2. Select and implement image 
metrics 


g. Process imagery and assemble 
database of metrics and 
corresponding ground-truth data 


4, Perform statistical analyses to: 
a. Identify ground-truth variables 
that most affect metrics values 
Develop method to predict 
metrics values given ground- 
truth data 


b. 


Figure 1. Overall methodology 


Table 1 
Sampling Locations and Times 


| Terrain Type Dats Collection Bite (6) 


Aberdeen Proving Ground, MD (APG) | 16-21 Mar 91 
18-16 Jul 4 


24-27 Mar 01 
10-13 Jul 91 


Fort Drum Military Reservation, NY 
(FTD) 


Fort A.P. Hill, VA (APH) 


15-18 Mar 91 
10-21 Jul 01 


9-13 Aug 91 
10-15 Sep 90 
Image data were obtained in visible, mid-, and far-infrared wavebands. 
Imaging systems were all off-the-shelf commercial systems; specifications are 
listed in Table 2. Thermal cameras are DC-restored calibrated systems. The 
visible camera is uncalibrated but uses fixed manual settings (i.e. no automatic 
gain control); therefore, comparisons can be made between images. Table 3 is 


an inventory of imagery that passed quality assurance tests and that were suc- 
cessfully calibrated as described by Berry, Rivera, and Sabol (1993). 


Coastal Plain Eglin Air Force Base, FL (EGL) 


Yuma Proving Ground, AZ (YPG) 
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Table 2 
Imaging Equipment Specifications 


Photometrics 200 
camera, with 


Thomson CCD 
datector 


Specification 


Agema Thermovision 
870 ayatem, infrared 
camera 


Agema Thermovision 
782 system, infrared 
camera 


Wavelength band 2-6.6y 
(no filter) 
3.8-5.6p 


(with SRX filter) 


FOV lens 2.5° x 2.8° 36° x 3.6° 3.9° (horizontal) 
x 2.6° (vertical) 
140 x 140 pixels 140 x 140 pixals 384 (h) x 576 (v) 
pixels 


Radiometric 8-bit resolution 6-bit resulution 14-bit resolution 
resolution 

Radiometric 0.1 °S at 3 °C 0.1 °C at 30 °C Not calibrated 
sanaitivity object temperature object temperature 

Radiomatric Not calibrated 
accuracy 


| Table 3 
Inventory of Imagery That Passed Quality Assurance and Calibration 
al 


= 
ra | Wavabond ine 
Ecco 


148/154 | 3077068 | 145/308 


173/178 | 301/208 | 106/308 


84/84 180/182 | 174/162 


112/112 184/196 134/196 


_ 64/66 153/154 | 0/164 
77/68 149/164 | 64/154 
FTD/7 

fea 


127/140 
Total 4,076 
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Weather conditions including air temperature, relative humidity, solar radia- 
tion, wind speed and direction, and precipitation rates were measured and 
recorded using a portable weather station interfaced with a digital recording 
unit. Additionally, radiometric temperatures of representative samples of the 
predominant terrain surface types were measured with staring radiometers. 
Automated measurements were made once per minute and stored as 15-min 
averages, Table 4 contains a description of the instrumentation used and the 
units of measurement. 


Terrain attribute data were measured for each scene through mensuration of 
color photographs taken from the designated sensor positions. Color 
photographs were taken of each 2.5° x 2.5° scene once during daylight hours 
using a 35-mm camera equipped with a zoom lens. Enlarged photographs 
were interpreted to delineate the boundary of each terrain/vegetation type 
within each scene. Boundary overlays were mensurated to determine: (a) the 
percent cover of each terrain type in the image plane and (b) the length of 
edge between all different types of terrain cover. The status of the vegetation 
(GREEN, active or dormant) was also recorded. 


In addition to the terrain variables determined by direct mensuration, 
several "generic" terrain variables were introduced to generate 2 terrain data 
subset common to all sites. Such variables include POLYGONS, TYPES, 
VEGET, EDGE, and HARDEDGE. POLYGONS consist of the total number 
of discrete enclosed terrain areas within a scene without regard to duplication 
of types. TYPES is the total number of terrain classes (Table 4, Report 1). 
VEGET consists of the sum of the percentage of grass cover (GRASS) and the 
percentage of tree cover (TREES). EDGE is the total linear length of edges 
(in units of angular degrees) between adjoining terrain classes. HARDEDGE 
consists only of linear edge between terrain surfaces with significantly different 
thermal properties that would be expected to have a strong thermal contrast 
during most times of the day. These include edges between water and all 
other surfaces, those between sky and all other surfaces, and those containing a 
vertical discontinuity such as trees and other terrain classes. 


Scenario variables were computed for each individual image, including 
camera pointing angles (AZIMUTH and elevation (ELEV)), range to the center 
of the image (RANGE), and bidirectional angle (angle between viewing and 
solar illumination vectors). Range to the center of each image was computed 
using a passive ranging technique applied to digital terrain elevation for each 
Site. These procedures are detailed in ECTA Report 1. Bidirectional angle 
(BDA) was computed based on the camera pointing angle and the solar illumi- 
nation angle at the time each image was collected. Each BDA data point was 
assigned to one of five classes: class "0" when the sun was not up, and 
classes "1" through "4" for angles between 0 and 180 deg in 45-deg 
increments. 
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Table 4 
Meteorological and Radiometric Instrumentation 


Datebase 
Variable Name* 
Fenwal electronic 20.4 °C Degrees C | AIR_TEMP 
UUT51J1 thermistor 
Phys-Chem. research REL_HUM 
PCRC-11 
humidity guage 
ft Incoming solar Li-Cor L1200S Watte/ SOL_RAD 
radiation (0.5- Pyranometer 
1.0) 
Met One, Inc., model WIND_SPD 
014a contact 
anemometer 


Met One, inc., model WIND/_DIR 
024a potentiometer 

windvane 

Texas Electronics, model | +1% at $2 

TES26 tipping bucket intr 

rain gage 


Everest Interscience, 

Inc., series 4000, 8 12u 

infrared temperature 

transducers 

(4"FOV) GRASS, ROAD, 
SOIL, TREE, 
WATER) 


All temporal (meteorological and radiometric variables), scene characteris- 
tic, and scenario data associated with each individual image were placed into a 
database for use in analyses. Table 5 contains a comprehensive list of these 
variables. 


Image Metrics 


Image metrics refer to the process and results of quantifying the distribution 
of specific features within a digital image. Image metrics are usually classified 
as to whether they require knowledge of target location (target dependence or 
independence) and whether they measure features in the entire image (global) 
or only a localized area (local). Numerous image metrics have been proposed 
and described (Beard, Clark, and Veltin 1985, ERIM 1985, Peters 1988). 

They range from very general global measures, such as image gray level mean 
and standard deviation, to very specific local measures, such as target-sized 
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Methodology 


contrast at specific image locations. Most image-based ATRs use contrast and 
edge strength of target-sized image features to identify regions of interest 
within the image (Peters 1988; Lahart, Jones, and Sheilds 1988). Research has 
shown some statistica! correlation between image metrics based on these types 
of features and detection/false alarm performance of ATRs. Because almost all 
ECTA imagery are without targets, only global target-independent metrics, or 
"scene metrics," are used. The following paragraphs describe these metrics 
and their computation. 


Twenty-one different scene metrics were computed in this study. These are 
listed and briefly described in Table 6; computational procedures and inter- 
pretation are described below. Twelve of these characterize the data-space 
distribution of temperature or brightness values; the remaining nine describe 
target-sized spatial variations within the imagery. Several metrics are 
dimensionless ENTROPY, SKEWNESS, KURTOSIS, REYNOLDS); all 
others have the units of the original image, which in the case of thermal 
images are converted to equivalent temperature units. The field of view (FOV) 
of the ground-obtained ECTA imagery often contained more foreground 
(objects at a closer range than 200 m) and far background (objects at a range 
in excess of 10 km) than would be expected with helicopter-obtained imagery; 
steps were taken to avoid these areas during image metrics processing, All 
metrics are computed only within a processing region of the image. This is 
defined as the central portion of the image remaining after the far background 
and foreground portions have been excluded. Passive ranging (Berry, Rivera, 
and Sabol 1993) was used to determine these areas. 


Data space distribution metrics provide general information about the statis- 
tical spread of digital brightness values but provide no information about how 
brightness values ale spatially distributed within the image. They were 
selected because of their common use as statistical descriptors. By and large, 
they provide no information about the distribution of targetlike features in the 
imagery. Processing was initiated by generating histograms of gray level 
values within the processing region. Eight measures are computed from the 
histogram: the minimum value (MIN), the 5-percentile value (PERC_05), the 
median value (MEDIAN), the mode (MODE), the 95-percentile value 
(PERC_95), the maximum value (MAX), the difference between the 95- and 
5-percentile (RNG_90), and entropy (ENTROPY). ENTROPY is a dimen- 
sionless measure which indicates uhe evenness of the distribution. Values near 
zero indicate an uneven distribution in which most of the pixels fall into rela- 
tively few different brightness values. High values indicate that pixels are 
evenly distributed among the available brightness bins. Since this measure is 
sensitive to variable camera sensitivity settings, entropy for thermal images 
was computed on 0.1 °C bins instead of raw digital values. The equation used 
is as follows: 
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i 
ENTROPY = ->> p, in (p)) (1) 


ie} 


where 
7 = number of brightness value bins 
p, = probability of a pixel occurring in the / bin 


The four moments of the distribution of gray level values within the pro- 
cessing region (mean, variance [standard deviation], skewness, and kurtosis) 
are computed (Press et al. 1986). The first and second moments, mean and 
standard deviation, have the dimensions of the original data—gray levels or 
equivalent blackbody temperature. Skewness and kurtosis are dimensionless 
numbers. Each is described below: 


N 
MEAN = © Y) 4, (2) 


a= 


N 
sp = |—1__ 3 (x, - MEANY @) 
N -. 1 ne) 


where 

N= number of pixels in processing region 

x, = brightness value of n pixel in processing region 
Standard deviation represents the root-mean-square deviation. 

Skewness, the third moment, is a dimensionless value which characterizes 
the asymmetry of the distribution. Values near zero represent a symmetric 
distribution. Positive values indicate an asymmetric distribution distorted to 


the right. Negative values indicate an asymmetric distribution distorted to the 
left. 
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3 
» |x, - MEAN 
SKEWNESS = 1 5x |2s 7 @) 
N 4 [SD 


Kurtosis, the fourth moment, is a dimensionless value indicating the relative 
peakedness or flatness of the distribution relative to a normal (Gaussian) distri- 
bution. Values near zero indicate a normal distribution. Positive values indi- 
cate a highly peaked monomodal distribution. Negative values indicate a flat 
or possibly multimodal distribution. 


4 
1 & |x, - MEAN (5) 
KURTOSIS =|. ees, | | 9 
URTOSIS | = | 


The remaining metrics characterize target-sized spatial variation in imagery. 
These were selected because they quantify the distribution of targetlike features 
and therefore may be relevant to ATR sensor performance. Two metrics char- 
acterize local variability within the processing region. Hetzler et al. (1987) 
described what they called a clutter metric (CLUTTER), computed as the aver- 
age standard deviation of adjoining squares twice the size of the longest 
dimension on the intended target. This was implemented using a square image 
window equivalent to 20 m on a side computed using the passive range esti- 
mation at the center of the processing region. This measure represents the 
average local variation. 


CLUTTER = EE Ee - 1.) (6) 
Moe OO” 


M = number of adjoining square windows in processing region 


where: 


J = number of pixels in square window 
Xmj = Stay level of the * pixel in the m’* window 
#, = Mean value of the m” window 


Reynolds (1990) extended this concept by computing the ratio of global to 
local variation. 
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SD - CLUTTER (7 


‘NOLDS = 
REY 3D 


REYNOLDS is a dimensionless ratio ranging between zero and one. Values 
near zero indicate that total variation (SD) is attributable entirely to variation in 
local regions (CLUTTER); conversely, values near unity indicate that local 
variation contributes nothing toward total variation. 


The remaining metrics characterize the distribution of target-sized local 
contrast. The basis for these metrics is computation of local target-sized con- 
trast using a double window (Carlson and Radford 1986). For each pixel 
within the processing region, a smallest enclosing rectangle is computed in 
image dimensions (pixels) corresponding to a 6-m (horizontal) by 3-m (verti- 
cal) rectangle! in scene dimensions (Figure 2). This computation is based on 
the passive range estimation of the current pixel. The rectangle, referred to as 
the target window, is surrounded by another twice its size, referred to as the 
background window. The background window surrounds the target window but 
does not include it. The mean brightness value of pixels within the background 
window is subtracted from the mean within the target window: 


tr D 
contrast = r XY x- 3 ¥x, (8) 


where 
T = number of pixels in target window 
x, = gray level of the ¢” pixel in the target window 
B = number of pixels in the vackground window 
x, = gray level of the 5” pixel in the background window 


The resulting local contrast value is stored in a contrast histogram. Positive 
values indicate bright (hot) contrast of the target window; negative values 
indicate dark (cold) contrast. This double window is convolved with all pixels 
in the processing region, resulting in a histogram of target-sized contrasts. 
Seven points are sampled on this histogram: the minimum (CNT_MIN), the 
5-percentile value (CNT_0S5), the 25-percentile value (CNT_25), the median 
value (CNT_50), the 75-percentile value (CNT_75), the 95-percentile value 
(CNT_95), and the maximum value (CNT_MAX). 


' These dimensions correspond to the smallest enclosing rectangle required to contain any 
armored vehicle from any forward-looking perspective. 
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ee 


TARGET WINDOW 


BACKGROUND WINDOW 


Figure 2, Range-dependent double window for measurement of target-sized 
contrast 


Metrics described above are listed and explained in Table 6. The image 
processing software for computing these metrics was written in Turbo Pascal 
5.0 and run on an IBM compatible 386-based personal computer equipped with 
a Paradise VGA Graphics card for image display. A listing of these programs 
is contained in Appendix A. 


Analysis Procedures 


Analysis was performed in two phases. In the first phase, graphical 
techniques were used to examine the range and variation of the ground-truth 
variables and image metrics. In the second, analyses were performed to deter- 
mine the importance of each ground-truth variable relative to selected image 
metrics and to use these ground-truth measurements as independent variables 
to predict metrics values (dependent variables), Several ground-truth variables 
thought to be important represent a nominal data type! (classes of bidirectional 
angle [BDA] and classes of plant growth status [GREEN)]); these data types 
cannot be directly used in multilinear regressions or other standard parametric 
statistical techniques. Consequently, an altemative non-parametric analysis 
procedure, "decision tree classifier methodology,” was selected. A decision 
tree classifier is a tree data structure’ which is traversed to arrive at a predic- 
tion of the dependent variable. The path selected at each node in the tree is 
determined by decisions on the value or class of independent variables. There 


1 Observations may be separated according to categories. 


2 Conventional data structure terminology, such as that found in Aho, Hopcrofi, and Uliman 


(1983), is used in describing tree str-ctures in this report, 
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Table 6 
image Metrics List and Description 


Minimum (MIN) & percentile 
(PERC_05) median (MEDIAN) 
mode (MODE) 95 percentile 

(PERC_(%) maximum (MAX) 


0% range (RNG_0O) 


Histogram-based measures of distribution of gray 
levels in image, units of °C in thermal imagery, units of 
digital brightness in visible imagery 


Difference between 95- and &-percentile values, range 
encompassing 90 percent of image, units of °C in 
thermal imagery, units of digital brightness in visible 
imagery 


Dimensionless measure of uniformity of distribution of 
pixel brightness 

Interpretation: high value (>3) indicates relative 
uniformity 


Entropy (ENTROPY) Carlson and Radford 


(1986) 


Mean brightness (MEAN) Average brightness in processing region,” units of °C 
in thermal imagery, units of digital brightness in visible 


imagery 


Standard deviation of all pixels in processing region, 
units of °C in thermal imagery, units of digital bright- 
ness in visible imagery 


Press ot al. (1986) 


Standard deviation (SD) Press at al. (1986) 


Skewness (SKEWNESS) Dimensionless measure of asymmetry of distribution of 
pixel brightness in processing region. 

Interpretation: negative values indicte asymmetric tail 
In negative direction; near zoro indicates symmetry; 
positive values indicate asyrnmatric tail in positive 
direction 


Press et al. (1986) 


Kurtosie (KURTOSIS) Dimensionless measure of peakedness of distribution 
relative to normal distribution. 

Interpretation: negative values indicate flat or 
multimodal distributions, vaiues near zero indicate a 
normal distribution, positive values indicate a highly 


peaked monomodal distribution. 


Average etardard deviation of boxes in image twice 
the longest dimension of target 

Interpretation: high values indicate local variation of 
pixel brightness in image 


Press et al. (1986) 


Georgia Tech clutter metric 
(CLUTTER) 


Hetzler et al. (1987) 


Reynolds metric (REYNOLDS) | Dimensionless ratio (0-1) equal to 
[(8D - CLUTTERVSD] 
Interpretation: represents the portion of thermal varia- 


tion attributed to variations within local regions 


Reynolds (1990) 


Target-sized contrast A range-dependent double-window matric which histo- | Sabol and Hall (1980) 


(CNT_nn) grams local target-sized contrast for all of image in 
nn=MIN processing region. n percentile represents the contrast 
08 brightness at the n percentile point on the histogram; 
26 x, 95 percentile contrast represents contrast which is 
80 exceeded by only 5% of image. 

78 Interpretation: high values indicate background con- 
05 tains target-sized features. 


MAX 


* Only pixels between 0.2- and 10.0-km range are processed. 
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are numerous types of decision tree classifiers; these are grouped based on 
type of tree structure, ways the tree is traversed, and methods by which the 
tree is "grown." Detailed reviews of decision tree classifiers may be found in 
several technical papers (Safavain and Landgreke 1991; Dattatreya and Kancl 
1985). For the present study, a hierarchical top-down binary decision tree 
methodology contained in the commercial software package CART (Classifi- 
cation and Regression Tree) (Breiman et al. 1984), was selected. Using a deci- 
sion tree analysis procedure made it possible to predict classes of a dependent 
variable (image metric values in this case) using independent variables (envi- 
ronmental measurements in this case) of unrestricted data types (nominal, ordi- 
nal', or interval’) without restrictive assumptions about the distribution of 
these variables. It further provides insight into how important each candidate 
independent variable is to the prediction. 


CART constructs a tree by recursively searching over all independent vari- 
ables to produce a sequence of optimal binary splits within the data set. A 
large tree is grown in this manner, subsequently, a pruning algorithm is 
applied which removes all branches within the tree which reduce overall accu- 
racy. The resulting tree maximizes predictive accuracy relative to its size. 
Figure 3 illustrates a hypothetical CART tree which predicts the class of a 
dependent variable. 


Te tree is entered at the top (root) node. Here a yes/no question is 
answered about the value of independent variable x(3). If x(3) is less than 5.0, 
the left descending branch is taken; if it is not, the right descending branch is 
taken. At each intermediate (nonterminal) node another yes/no question is 
answered regarding the value of other independent variables (x(n)), dictating 
the path to be followed. Arrival at a terminal node ends the traversal process 
with a prediction of the dependent variable’s class. The tree is grown and 
tested in a two-step process, beginning with the random division of the data set 
into two portions. One portion is used to grow the tree. The accuracy of this 
tree is then determined by processing the second portion through the tree and 
comparing predictions of the dependent variable against corresponding actual 
values. In this study the dependent variables are binned into ordinal classes so 
predictive accuracies are computed from the resulting confusion matrices. 


Confusion matrices, also known as error matrices, tally the predicted class 
of each individual observation against its corresponding actual class in a tabu- 
lar form. The number at each location in the table (row i, column j) represents 
the number of observations which were classified as class i but were actually 
class j. Table entries along the main diagonal (i=j) represent correctly classi- 
fied observations; those off the main diagonal represent specific misclassifica- 
tions. "Classification accuracy" is computed as the sum of the main diagonal of 
the confusion matrix divided by the grand total of the entire matrix. Since the 
dependent variables are of ordinal data type, all misclassifications are not 


1 Observations may be arranged from smallest to largest. 
2 The numerical value of the observation has physical meaning. 
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Figure 3. Typical binary tree structure generated by CART 


weighted equally. Misclassification into a neighboring class (e.g. classifying a 
true class #1 as a class #2) is a less severe error than classifying into a no 1- 
neighboring class (classifying the same true class #1 as a class #3). A second 
accuracy measure, "severe misclassification," represents the percentage of the 
observations misclassified into a non-neighboring class. 


CART internally evaluates the relative importance of each independent 
variable by quantifying how well data splits on each independent variable 
separate the dependent variable into distinct groups at each node within the 
tree. This measure is summed over all nodes for each variable and normalized 
by dividing by the highest value produced among all variables. The result is a 
range of importance values between 0 and 100, with the most important vari- 
able always having a value of 100. This procedure, described in greater detail 
in Breiman et al. (1984), is used to judge how important each independent 
variable is toward the prediction of the dependent variable. 


Several of the metrics computed quantify image attributes similar to those 
used by ATR systems in th !»cation of regions of interest. These include 
Standard deviation (SD), CLUTTER, REYNOLDS, and the local contrast 
metrics (CNT_nn). Of these, SD and CNT_95 have demonstrated a statis- 
tically significant correlation with probability of detection tor several actual 


Chapter 2. Methodology 17 


18 


ATR systems’. Similar data do not currently exist to support the importance 
of CLUTTER and REYNOLDS. Relatively high values of these selected 
metrics tended to coincide with relatively low values of probability of detec- 
tion. For the purposes of this study, ordinal classes of SD and CNT_95 in the 
visible and the 8-12» wavebands” serve as the dependent variables for which 
predictors will be developed. Ordinal classes were derived by histogramming 
these metrics and binning each into low (0-33 percentile), medium (34- 

66 percentile), high (67-89 percentile), and very high (90-100 percentile) 
classes (Figure 4). These ranges are considered representative of values for 
these sites, and binned classes are used as indicators of scene complexity. 
CART trees were generated to predict these metrics classes. The question of 
importance of the various ground-iruth variables is directly evaluated by 
CART. 


Decision tree classifiers were generated for classes of metrics SD and 
CNT_9S in visible and 8-12» wavebands using three separate sets of indepen- 
dent variables as illustrated in Table 7. The first set (level 01) used the full set 
of scene, meteorological, and scenario variables. The second set (level 02) 
used a reduced set of generic scene data, full meteorological data, and scenario 
data. The third set (level 03) used only meteorological and scenario data. 
Meteorological data sets used were the same for thermal and visible bands 
with the exception that only instantaneous solar measurements (as opposed to 
time series measurements) were used for visible analyses. The relative impor- 
tance of independent variables arid the overall accuracy of resulting classifica- 
tions were compared to determine the importance of the various data types. 


1 This information is contained in the final briefing of the MultiSensor Fusion Demonstration, 


presented at Fort Eustis, VA, August 1989; available through DTIC. 
2 Data losses in the 3-5p waveband imagery (Table 3) preclude extensive analysis in this band. 
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Figure 4. Four classes of dependent variables: (a) thermai SD, (6) thermal CNT_95, 
(c) visible SD, (d) visible CNT_95 
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Table 7 
Classification Trees Generated for Metrics SD and CNT_95 for Visible and 8-12 
Waveband imagery 


Dependent Variabie 03) Meteorological No 

Groupings (Waveband) | Scene Scenario Generic Scene Scenario | Scene Data Scenario 
[Therma tauos.cratios | stauioz,ctauoz | stautos, cTauios | 
[Visbie dL svatios, cvaior | svatioz,cvaiioe | svalios, cvaLtos | 


TREE FILE LEGEND: mwsssw - 7 character code 

where 

m: metric selected (8 = SD, C = CNT_O5 

w: wavelength (T = thermal, V = visible) 

ss: site (ALL = all sites) 

wy: independent variable data set 

(01 = meteorological, scenario, and terrain data 

02 = all meteorological and scenario data with generic torrain data 
03 «= meteorological and scenario data only) 
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3 Results 


Comparison of Ground-Truth Conditions and 
image Metrics 


Scene and scenario data 


A complete listing of data from the scene content analysis, mentioned in 
Chapter 2 and described in detail in Report 1, is contained in Appendix B. 
Major season independent attributes of the scenes are illustrated in Figures 5 - 
7. Terrain classes comprising each scene at each site are shown in Figure 5. 
Number of discrete terrain areas (POLYGONS) for each: scene at each site is 
illustrated in Figure 6. EDGE and HARDEDGE for each scene at each site are 
shown in Figure 7. 


Eglin AFB, Fort Drum, and Fort A.P. Hill sites contained similar types of 
vegetation cover, All were predominantly vegetated with grassy fields in the 
foreground and trees in the background. These three sites also had a relatively 
small portion of IIARDEDGE. The Aberdeen Proving Ground site was the 
only one to contain surface water, resulting in a relatively large portion of 
HARDEDGE. Further, the Aberdeen site had the largest number of discrete 
terrain areas (POLYGONS) among all sites. The Yuma Proving Ground site 
differed most from the other four, it contained the largest portion of bare 
ground and nonvegetated terrain and was the only site with large surface 
geometry variations (i.e., mountains). 


Meteorological and radiometric data 


A complete listing of meteorological and radiometric data is contained in 
Appendix C. Temporal plots of meteorological and radiometric conditions 
during the cool-weather excursions are illustrated in Figure 8; those for warm 
weather excursions are illustrated in Figure 9. During the cool weather 
excursions (Figure 8), cloudy conditions existed at the Fort Drum and Fort 
A.P. Hill sites; clear sky conditions were present only at the Aberdeen site. 
Accordingly, air temperature was warmer and exhibited greater variations, 
relative humidity was lower, and thermal contrasts between different terrain 
surfaces were higher at the Aberdeen site. During the warm-weather 
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Figure 5 
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Figure 6. 
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Figure 7. Total linear edges (EDGE - total height of bar) and edges between thermally dis- 
similar terrain classes (HARDEDGE - thickness of darkened portion of bar) by 
scene for each site 
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Figure 8. Temporally varying meteorological and radiometric conditions by site for cool-weather excursions (Continued) 
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Figure 8. (Concluded) 
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times and conditions. Table 11 lists the classed outputs from this algorithm 
for these conditions. Predicted class levels for the thermal local variability 
metric indicate low levels from midnight (0000) unti] 10 AM (1000). Begin- 
ning at the noon (1200) measurement, the levels jump to high then fluctuate 
between medium and high until 10 PM (2200). With these predictions as a 
guide, critical test times can be selected to achieve the best balance of scene 
complexity conditions available at this site and time. The predictor could also 
be used in a site selection mode by entering characteristics of potential scenes 
and expected seasonal meteorological conditions to determine which candidate 
scenes give the best range of complexity conditions for the intended system 
mission and operational environment. 


As with any empirically based research effort, the limits to which the 
results may be safely generalized are based on the adequacy of the sample and 
the underlying assumptions. Practical considerations dictated that the sample 
from which analyses were performed be of a limited size and be obtained by 
systematic sampling as opposed to random sampling. All imagery data were 
obtained from a single sensor location at each of five separate test locations. 
These were some of the most commonly used ATR development and testing 
sites; however, they represent only a limited sample of terrain conditions. 
These included three temperate locations with deciduous vegetation, each 
visited once during leaf-on and leaf-off periods, and one coastal-plain and one 
semi-arid location, each visited during the summer. Temporally varying 
meteorological measurements also represent a limited range of conditions. The 
systematic diumal sampling was intended to capture all conditions occurring 
within the 24-hr period sampled, but the 24-hr period is essentially just a sin- 
gle random sample of all daily weather conditions occurring at a site during 
the season sampled. Completely sunny conditions were encountered during the 
Yuma visit and during both Aberdeen visits. Likewise, completely cloudy 
conditions were encountered at Fort Drum, No snow conditions were 
encountered at any site, Clearly, there are myriads of important meteorological 
conditions and terrain factors not considered in this study because of the 
limited database. It is doubtful whether the present database could be consid- 
ered an adequate sample of all conditions and terrain factors likely to be 
encountered at CONUS test sites. 


The second consideration in generalizing the results is the degree to which 
the selected metrics really reflect difficulty posed to an ATR. Previous tests! 
showed a high negative correlation between probability of detection (averaged 
over all targets in a scene) and CNT_95 and SD fot some or all ATR systems. 
CNT_95 showed the strongest correlation (p<<0.01) with average probability 
of detection for all ATR systems in that test”, However, even 


1 Information contained in final briefing of the MultiSensor Fusion Demonstration, presented at 
Fort Eustis, VA, August 1989; available through DTIC. 

2 This strong negative correlation indicates that high values of CNT_95 terd to coincide with 
low values of probability of detection averaged over all targets in the scene. 
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Table 11 
Predictions of Local Thermal Variability Using Decision Tree 
CCTALLO4 


Predicted Level of Local Thermal 
Time of ae (hour) — 


this measure explained only 5 to 17 percent of the overall variation in average 
probability of detection. Factors such as target type, target orientation, opera- 
tional condition, and specific location of a target within the scene must be con- 
sidered, The predicted class of the metrics should be taken only as a 
"broadbrush" indicator of conditions. When indicator metrics values are rela- 
tively low, indicating a bland scene, probability of detection will on the aver- 
age tend to be higher than when indicator metrics are relatively high. 


Given these considerations, what uses and limitations are recommended for 
the resulting capability? As with any empirical predictive technique, predic- 
tions are generally satisfactory if inputs are within the range of values used to 
develop the predictor. That is, the predictor should be used to interpolate not 
extrapolate. Metric class predictions for additional scenes from the five imag- 
ing locations under similar weather conditions would expectably be as accurate 
as those of the test data sets (Table 9). Predictions for conditions outside of 
the range of those in the database may have lesser accuracy. 


Taking an empirical approach toward solving this complex problem was 
motivated by the belief that a first-principles approach would have involved 
attacking overwhelming levels of complexity (Sabol and Hall 1990) with the 
limited resources of this study. The robustness of the resulting empirical pre- 
dictive technique is unknown but is suspect based on the limited database size 
consideration described above. The standard solution to this problem is to 
increase sample size. However, we feel that this approach would yield only 
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limited return. Field data collection is costly and labor-intensive, and the prob- 
lem may be too complex to adequately solve with a strictly empirical method. 
Efforts are currently underway within the Department of Defense to better 
understand this problem from a first-principles approach. The tri-service Smart 
Weapons Operabiiity Enhancement Program (SWOE) managed through the 
U.S. Army Engineer Cold Regions Research and Engineering Laboratory, 
Hanover, NH, and supported by WES and other agencies, is developing first- 
principle computerized scene generators for passive thermal and active milli- 
meter wave (MMW) sensor systems. Validated thermal and MMW scene 
models are scheduled for release by fiscal year 1995. Coupling such a model 
with a sensor transfonnation model would produce a synthetic image which 
could then be processed through an ATR logic set to simulate system #erfor- 
mance directly. Sensitivity studies using such a capability would enavie the 
analyst to evaluate the effects of selected scene formation factors directly. 
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Al 


PROGRAM: THERMAL.PAS 
MAIN PROGRAM FOR 
THERMAL IMAGES 


AGAUTILS,PAS 


CALIBRATION 
EQUATIONS 


EAD FIRST 
SITE TO BE 
PROCESSED 


READ 5 
MAGE NAME TO 
BE PROCESSED 


READ NEXT 
SITE TO BE 


ROVESSED 


READ IMAGE'S 
AZIMUTH ANO 


AGARANGE.PAS 


AGAUTILS.PAS 
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(i 5 cam tne 


AGAUTILS.PAS 
AGACOMPT.PAS 
AGAMETH.PAS 


IMAGE 
METRICS 


| COMPUTE MORE 


DISCARD 
IMAGE 


CONVERT IMAGE 
METRICS FROM . 


PRINT OUT 
IMAGE 
M 
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A3 


~~ 


(tue) 
Program Thermal; 


ees 

ert, 
ageutila, 
aqararge; 
a, 
printer, 
egecenpt, 
agemetre; 


omalligly 
large_ely 
teep ain 
temp_osk 
temp_nod? 
temp_05 
tenp_95 
teap_med 
temp_meen 
temp_ed 
eeen 
entrapy 
etd_dev 
okewnean 
hurteshs 
upbound 
\exbound 
GLW_hiat 
luaet ile 
firat_caleq 
lagtile 
outtile 
tap_neme 
data_loc 
ent 

code 

view img 
view 

axel file 
ardeg 
stain 


iuage_name; 


integer; 
GLV_hlato; 
text; 
caleqptr; 
text; 
Rent} 
ating (1133 
atringi2; 
integer; 
integer; 
integer; 
integer; 
text; 
integer; 
integer; 


THERMAL 1 OF 4 


A4 
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azsec 2 integer; 
eldeg 1 Integer; 
elmin t Integer; 
elaec + integer; 
e i POLWTER; 
slres_file t tent; 
{ege_filename : string(t2); 
matr_tilenmme : steing(!2); 
etre : stringl2); 
windows + window_erray; 
clutter s reel; 
feynold t real; 
histo t Bisteptr; 

: 

t 

t 

t 

, 

t 

t 

t 

t 


THERMAL 2 OF 4 


cont_ing cont_arrey; 
wincont integer; 
maxcont integer; 
theent_05 real; 
theent_25 real; 
theent_$0 reals 
theont_75 reat; 
theont_95 real; 
leg integer; 
0 ergrinteger: 
usbyte;¢) 


begin 
READ_CALEGR( Firat celed); 


apsign(eites_file, ‘SITE_ALL.IXT!); 
reset(aites file); 


readin(sites file, imgs_tileneme, atr2, metr_filensme); 
while( (mgs_tilename «4 ) do 
begin 
ansign(iegfile, 16G8_FILENAME); 
roset(imgtitey; 
readinilmgtile, dete_loc); 
ent rm 0; 
resdin(imgtile, tap name); 
while( twp_name <>‘! ) do 
begin 
eseigniouttile, META_FILENAME); 
($1- 
eppend({outfiia); 
cate} 
44 (LOreeult <> 0) then rewritecoutfile); 


ent :# ent ¢ 1; 

feme 1* copy(tep_name,1,8) ¢ '.' + copy(tep_neme, 10,2); 
writelnceutput, "Processing image @', ents4, ' ', data_loc, name); 
GetDeted junk, procaonth, procday, junk); 

GetT imetprochour,procain, Junk, junk): 


{¢ copy(name,1,2) = 'EG' then 
easigniaze! file, "EG_AZEL. TXT!) 
elee 
{f copy(neme,1,2) = ‘YG! then 
aseignceael file,  *YG_AZEL.THT') 
else 
1f copy(neme,1,2) = ‘AH! then 
ecnign(azet_file, ‘AH_AZEL. TXT‘) 
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else : THERMAL 3 OF 4 
{¢ copy(name,',2) © "AG! then 
anelentazel_file, *AG_AZEL. TXT) 
else 
if € copy(name,1,2)> = ‘f0) and ¢ copy(nems,11,1) 2 ‘1! ) then 
aesiontazel_file, ‘FO_AZBLI.TXT!) 
else 
{¢ ( copytneme,1,2) © /FD') and ( copy(name,11,1) = '2! ) then 
aasiontazel file, ‘FD AZUG2, 1X1’) 
else 
{f copy(neme,1,2) = HL’ then 
espigntaxel_tile, ‘ML_AZEL.TXT*); 


IMT T WINDOWS (windows, name); 
Pevet(azel file); 
repeat 
val (copy(neme,7,2), viewing, code); 
poadiniazel file, view, azdeg, amin, atsec, eldeg, elmin, elsec) 
until ¢ view = view {mg >; 
closecarel_f (le); 
animuth t= axdeg > (((eztec / 60) + amin) / 60}; 
elevation 19 eldeg + ((delsec 7 60) ¢ alain) / 60); 


READ RANGE(date loc, nome); 


fo) (f ¢ Conpy(rame, 1,2) = YO") oF (copy(name,1,2) = ‘yg’) ) then 
FIND BOUNDS YPGCezimth, elevation, \evbound, ubound, name} 
else 
FINO BOUNDS ( lowbourd), Upbourdd, name); 


{f ((lowbound- upbound)<15) then 
begin 
UritELMC bounds restriction’); 
WALT! 
end; 


MARK(p)? 

READ_IMAOE (image, name, GLV_hiet, min, max, andr, perctS, perclS 
» Medien, mean, entropy, lewbound, Uupbound 
a firet_caleq, date_loc): 


Gmail gly se ( (gly mistcOysgiv miagiy+@ly ies (2) ely miet(296) > © 100; 
Carge_oly te ¢ Colv_nieti253)ogtv Niet (254}eqiv hiet (255) :/aiv histl2$6) > * 100; 


(6 ¢ (oma l_givetarge_giv) < 0.50 ) Then 

begin 
temp_ain tm GLV_TO_TE@ (ain, name, ¢iret_categ); 
Cemp_mex 1@ GLV_TO_TeWe(mex, name, first ealeg); 
temp_wode 1© GLV_TO_TEMP(mode, name, firet_caleq); 
temp_05 16 GLV_TO_1H@(perc0S, name, first caleq); 
temp_%S 1* GLV_TO_TEMPCperc®S, name, firet_caleq)) 
temp_med :© GLV_10_TEMP(medion, nem, firat_caleq); 
Cemp_mveni=s QLV_TO_TEM>Cround( mean), Name, first caleq); 


COMPUTE _STATS( (ange, etd_dev, min, perctS, mean, skewness 
« kurtosis, lowubound, upbound); 


temp_ed 1¢ DELTA _TEMP(round(std_dev), meme, firet_caleq); 


COMP _REYNLD( image, meme, upbound, loubound, clutter, reynold, atd_dev); 
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if clutter < (9999.9) then THERMAL 4 OF 4 
eiutter r= DELTA_TEMPCroundiclutter), name, first caleq); 


Wew{hlato); 


CONTRABT2(histo, lowtveund, uptourd, fange, windows, cont_img, nome); 
aimcont i= PERCENTILECKi ate, 1,0/1NT(histe*. bine (256) )); 
@axcont te PERCENTILEChiato, 0.9995); 
theent_O5 := OELYA_TERPCPERCENTILECHIsto,0.05), name, firet_caleq); 
theont_25 se DELTA_TENPCPERCENTILECHisto,0.25), name, firet_caieq); 
theent 30 10 OELTA_THNP(PERCENTILE(hieta,0.80), name, firet_caleq; 
theent_75 i= DELTA_TEMPCPERCENTILE(hieto, 0.75), name, first_eslea); 
theent_95 r= DELTA_TEMPCPRACENTILE(Hiato,0.95), name, firat_caleg); 
fleg rv 6; 
ond 
else 
begin 
teap_ain ie 99.9; 
temp_asn 1” 99,9; 
temp_mode 16 99.9; 
temp_05 1 09.9; 
temp VS 1s 99.9; 
tenp_aed t= 99.9; 
temp_msani= 99.9; 
temp_ed i= 99.9; 
elutter := 99.9; 
peynold 1* 99,9; 
theont_OS 1° 99.9; 
theont_25 15 99.9; 
theent_50 12 99.9) 
theont_75 1° 99.9; 
theont_93 1° 99.9; 


ord; 


flag te f 

ond; 

RELEARE(p); 

weitelntouttile, name:t0, ' ’, procdayié, ‘ ', precmonth:3, ' ‘ 
+ Ptochout:s,procminss, ' ', artmuth:10:4, ¢ ' 
1 Slevation: 10:4); 

writelntoutfile, upboundiS, § '. inwbound:5, ' ', temp_mgan: 10:3 
oft, templeinntOrd, 6 4, tempmanntors, ¢ ! 
» Comp_moder 10:3, ¢ *, omall_giviéi2, / 
+ Varge_glyi6i2); 

writeintoutfile, temp_05:10¢3, ¢ ', temp.05:1013, ' ', temp_medi t0i3 
+ ft, Remp_eds1013, ' 4, antropy:10:3, ¢ ' 
> akewnensi 1013, ° ', kurtoete:10:3); 

weitein(outtite, cluttar;1003, ' ', reynold: 10:3, ' *, theont_0511013 
a ff, *heane 2511013, ¢ %, theont_5011013, ¢ ! 
» Sheomt_75:10:5, / *, theamt 9511013, / ', flegity; 

Q 
elowecoutfiles; 


readiniimgfile, tap_name) 


ord) 
Closet imgfile); 
reedinieltes file, 


closecaites_ file) 


ond. 


fge_filename, sir2, metr_ filename} 
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Unit AQAUTILS; 
interface 


wees 
dos, 
printer, 
erty 


const 
hunrew © 140; 
fumcol = 140; 
left junk = 2; 
right. junk 2 3; 


type 
\mageptr © “imegerow; 
\maperow © record 
Gate : array(!..camcol) of byte) 
ond} 


{moge_fame u atring(12); 
videoptr = “videorgu; 
Videorew = reserd 

ond} 


centptr © “contrew; 
tentrew » recard 


ond; 


OLv hate 
Otringls © stringlt2); 


caleqptr © “ealeqe) 
oal_eqs § record 


ome t integer; 

\mg_name + etring(¥) 7 

unit_type oo: ateing(3); 

slepe 1 realty 

intercept + real; 

next _eateq | ealen ptr 
end) 


«* ww eeen Cs 
procedure read_imege(var image i 

name t 

ver Alete t 

war ain t 

ver max t 

var node i 

var perc 1 

ver percds 1 
1 

t 

t 

H 

1 

H 


ver median 

ver mean 

var entropy 
le bound 


upboursd 
tiret caleq 


As 


(nege_erray © array(1,.rumrew) ef imageptr; 


diop 1 array({..mumeol * 3) of syte; 


video array © arrey(l,. mare © 3) of videoprr) 


date: array(l..mumol) ef integer; 


eont array © array(1..ramrou) of contptr) 
® arrayl0..256) of integer; 
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ee | 
\mage_errey; 
image name; 
GLYN eto; 
Integer; 
integer; 
integer; 
integer; 
integer; 
imteger? 
real; 
real; 
Integer; 
integer, 
coleq_ptr; 
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Proccdure MEAD IMAGE INFOC name 1 
var day + 
var month | 
var year | 
ver hour ¢ 


var el 1 
Procedure READ BABEL INE INFOC nome 
var day | 
var month | 
vor year | 
var hour 
var minute: 
ver ext 
var el H 


(reer emwrerere racer eereennctaea tse etuaseane 


lop | erent ation 
ver 
range =| double; 
level 1 double; 
a 1 double; 
b t double; 
€ t double; 
1 otringis); 


Nrocedure read_image(var image 


var histo 
ver win 
var wan 
var anode 
ver pereds 
var percdS 
ver median 
ver mean 


dataloc _ 1 eteing 22; 


Procedure COMPUTE_STATS( (mage 1 lmage_arcay; 

war std_dev : real; 
ain t Integer; 
man t integer; 
mean 1 realy 

var skewneas 1 real; 

war kurtosis t real: 
\ewbound | integer; 
Upbound : integer); 


Proceture READ _CALEQS( var caleqi 1 caleq ptr); 


function GL¥_ to tenp( value 1 byte; 
jmagename 1 image name; 
flret_caleq | caleq ptr) i reat; 


function OELTALTENOCOLY t integer; 
imgneme t inage_neme; 
firatlealaq : caleq_ptr) 1 real; 


image, nemes 
integer; 
integer; 
inteper; 
integar; 
inteper; 
realy 
real); 

ti age name; 
integer; 
integer; 
Intwger; 
Integer; 
integer; 
real; 
real); 


see eenenweeey 


eosuaeeedeeeRe: 
{ tmage_erray; 
image_name) 

GLY histo; 

integer: 

integer; 

integer; 

integer; 

Integer; 

Integer; 

real; 
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AGAUTILS 


RUSUADEND POURSeR ECU Baten westudMennHE) 


AQ 


(evecrecmeccvoncsensenssaucons 


procedure HEADER(datal 


type 


Overflow 


A10 


var entropy : real; 
lowbound + integer; 
upbound t integer; 


flratcaleq + cateqptr} 
dataloc 1 atringt2); 


oo ot otringt2; 


ancovnonen) 


filename 1 imege name); 


apa_headersrecord 
Image_version 1 integer; 


integer; (eee manipulate) 
array{1,.é4) of byte; 


Propram_version t array(1..6) of char; 
origin t integer; 
lange id t arraylt..13) of char; 
dual 1 byte; 
tmoge, drive 1 erray(\.0a) of char; 
cuany 1 array(1..36) of byte; 
Carment t wrray(1..401) of wher; 
Prepresess + integer; 
Mag_tactor 1 double; 
Signal _eap t (ateger; 
Hanipul ated t Integer; (boolwan, but need to read 2 bytes) 
Qnapshat : Integer; Coee Kanipulete) 
Oi fterence t integer; (see marapulate) 
cum + byte; 
date : arrey{!..3) of integer) 
Time t arreytt..4) of integer; 
Title t array, 46) of chor; 
teamer 1 array(,.43) of char; 
and t oyte; 
Lene t arraytt,.43) of char, 
am 1 byte; 
a 1 double; 
ft] : double; 
' 1 double; 
Deal 1 double; 
ips t double; 
beta 1 double; 
dievel + double; 
tmiasivity 1 dovole; 
Ares_omias 1 double; 
Object_distance 1 double; 
Computed trans t double; 
Gatimated trans 1 double; 
At@oe temp 1 double; 
Ambient temp t double; 
DI tp 1 double; 
Olftq 1 double; 
a) t erray(1,.64] of byte; 
Level 1 double; 
fenge t double; 
Under f 1 ow 1 integer; (eee manipulate) 
! 
! 


1 file of age_tesder; 
| age_header; 
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AGAUTILS 


3 OF 10 


assign(infile, dataloc + filename); 
reset(infile); 
read(inftle,a header); 
range:se header .range; 
lavet:e0 header, level; 
a 1m » header.r 
b := e_header.b; 
¢ 38 e header. f; 
Scanner ¢= copy(e_heeder.scanner, 1,5); 
cloaecinfile; 
and; (header) 


(iseesee we eerererecererscensscssenecee 


function GLV_PERCENTILE( Ver histiGL¥_histo;percentrreal >): integer; 


var 
ttop_at : integer; 
{ t integer; 
aceum =: longint; 
» begin 


etop_at 2% round(percent®(mt(hist (256) )); 
eccum 38 0; 
{ 0; 
while (accum<stop_at) do 
begin 
eccum se accun> hist (i); 
frefeh 
end; 
OLV_PEACENTILE :® 4; 
ond; 


(cv vceveccens Peer errrrrrrrrre rere rire eecenenees Pere retry | 


var 
infile + file of byte; 
(,),% 1 integer; 
accum 3 longint: 
histed + arrayld,.512) of integer; 


prob: real; 
junk} byte; 
atart, 


current: integer; 
freqmax 3 integer; 


begin 


HEADER (dataloc, mama); 
assign(infile, dataloc + name); 
resetcinfile); 
for kse to 846 do 
reed(infile, junk); 
for k t# 1 to mumrow do 
mevt( image (k}); 


i:m2; (toed image) 

repeat (even rowa (lat interlace) 
for jr 1 to mumcal do read(infile, imageli}* detalii); 
ireie2; 

until (44 Crasirow #190; 

{uty 

repeat (odd rows (2nd interlace)? 
for jr: 1 to mumeol do read(intile, imaget!)*.dataly]}; 
(reledg; 
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AGAUTILS 


OF 


10 


All 


until (1 > rumrow); 
close(infile); 
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win 12 255; (compute image state on renge bounded image) 
@ax 18 0; 
eccun te 0; 
for | 16 @ to 286 do 
Aisto({) s= 0; 
for i:eupbound to Lowbound do 
for Jrecleft_junk+l) to (numeol-right_junk-1) de 
begin 
{¢ (image tii *.detat]) » mex) then 
max ve (mageli}*.date()); 
1f Cimeget{)*.detalj? < min) then 
min = inage(i)*.date(}); 
eccum :« sccum + imege(!)*.cete(}); 
hintolimegelt)*.datats?] := histoCimage({)*.datels}) + 1; 
hiato(2Sé) := histe (256) + 1; 
end; 
eon «18 accum/histol256}; 
pereOS 1° GLV_PHEACENTILE(histe,0.05); 
Median 1© GLV_PHRCENTILE(histo,0.50); 
perc9S 19 GLV_PERCENTILE(histe,0.95); 


atartse trunc(S*GLv_ To TEMO( min, name, first caleq)); entropy computation) 
tor ise to 512 do histozii} 190; 
freqmax := 0; 
for { :* min to mex do (compute mode and create thermal histogram with 0.2 deg C bin size) 
begin 
{f hiatol!) > treq_max then 
begin 
freqmax t= histolil; 
Mode 18 | 
ordi; 
current rm teune(S*GLv_TO_TEMPC | , name, firat_caleq)); 
histodtcurrent-stert) :# histoli) + histo2(current-start); 
ond; 
entropy :# 0.0; 
tor {18 0 to 512 do 
begin 
prod te hinto2(i)/histo{256); 
4f (preb>0) then 
entropy := entropy + (probd*in(prob)); 
end; 
entropy i* -entropy; 
end; (read_binary, image) 
(tucqeaeeu seuss eesuesnans cusses scenneene teenaeRe es ae se Reet eR EE anSaesauseEEs ) 
Procachire COMPUTE STATS( (mage + imepe_arrey; 
ver std. dev : real; 


ain t integer; 
on + integer; 
meen t real; 


var skewness : ceal; 
var kurtosis ; reel; 
lowbound : integer; 
upbound : integer); 


ver 


i, 1, & Los integer; 
eccuml s real; 
eccum2 : Peal; 


At2 


Appendix A Metrics Image Processing Sottware Source Code 


accund > real; 
square t real; 
temp 3 real; 
rampixel 1 real; 
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begin 
accual 12 0; 
accum2 18 0; 
eccumS 38 0; 
rampixel r= 0; 
for { 1= Upbound to lowbound do 
for | :@ 1 to mumecol do 
{f (Cj > heft junk) and ¢) < mumcol - right junk)) then 
begin 
rumpixel :© mnumpixel + 1; 
temp :0 image(i}*.datetl} - mean; 
sqpare 12 aqr( temp); 
aceunt i= accum! ¢ square; 
accum? 12 accum2 ¢ (squere * temp); 
eccunS 19 accumd + (aquere * square); 
end; 
accum) te accuml / (rumpixel = 1); 
std_dev :© sqrt(accua!l); 
eccum2 := accum? / (accuml * etd dev); 
skewness := accum2 / rumpixel; 
eccumS := accumS / (eqr(accum!)); 
kurtosts :=# (accumd / numixel) - 3; 
and; 
irri tie by 
Procedure READ CALEQS( var caleq! + caleqptr); 
var 
infile text, 
curr_cateq : caleq ptr; 
blank} string(3); 
blank35 string(35) 
blenk32 + string(52) 
obs : string(3); 
\mg_name atring(9); 
unit_type : string(3); 
slope : string(6); 
Intercept etring(7); 
code integer; 


begin 
new(calegt); 
curr caleq :* caleql; 
apsignCinfite, ‘cal_eqs. txt’); 
reset(intila); 
teadin(intile, blank}, obs, blank}, img_name, bi snk35 
, UNit_type, blenkS2, slope, blank3, intercept); 


while( abs <2 '' ) do 
begin 


val(obs, curr caleq*.obs, cade); 
curr_caleq™.img_name 12 img_name; 
curr_celeg*.unit_type :* unit_type; 
valislope, curr_caleq*.slope, code); 

val intercept, curr_caleq”. intercept, code); 
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weet es ete ee aang 


readintinfile, blenk3, obs, blank’, img_neme, blank3$ 
+ Mit type, blankS2, slope, blank3, {ntercept): 


if ¢ obe @ /4 > then 
begin 
new(eurr_caleq*.next_caleq); 
curr _celeq t= curr_caleq’.next_ealeq 
end 
alse 
curr_calaq*.mext_ealeq te WIL; 


end; 


closecinfile); 


end; 


(2neenseansasneeceniensenenseenessseeRnenns: RUCECOsee Re ceesneeeenesEsesERES) 


function GLV_to_Temp¢ value 1 byt 
imageneme 1 image name; 
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ver 
feetherne 
curr_caleg 
eq_feund 
obs 
unit_type 
slope 
intercept 
tone 
atring? 
inftle 
eg | nome 
egl_time 
egl_teap 
eo siv 
es\_rne 
6g! slope 


firet_celeq : caleqptr) : real; 


real; 
caleq ptr; 
boolean; 
integer; 
atring (3); 
real; 
real; 
real; 
atring(9); 


£ text; 


string (11); 
real; 
real; 


* resol; 


real; 
recl; 


begin 
IF ¢ COPY(imegername,1,2) «2 ‘BG’ ) then 
begin 
{aotherms :© level + range*(127-valwe)/256; 
(f isotheras < 1.0 then (sotherms := 1.0; 
curr_celeq := firet caleq; 
repeat 


Btring? 1° copy( imegename,1,6) * copy. imagename,9, 11); 
eq_found := ( curr_caleq*.iag_neme © string? ): 

ome 1= curr_caleq*.obs; 

un(t type te curr coleg*.unit_type; 

slope r= cure_caleq*. slope; 

Intercept r= cure _caleq™. intercept; 

curr_ealeq := curr celeq* neat caleq 


until( (curr_celeq © WIL) of eq found ); 


if not(eq_ found) then 
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begin 
writeln; 


hatt 
ond 
else 
begin 


else 
begin 


end; 
end; 
end 
else 
begin 


egl_neme re! fy 
reset(infile); 


closecinfile); 


end 
end; 


function DELTA_TEMPCGLY 


begin 
1 (GLV<0) 
then OELTALTEMPre © (QLV_TO_TEMP(abs(GLV), fegname, first caleg) + CLV_TO_TEMP(O, imgname, first_caleq)) 

@lee DELTA_TEMP:= GLV_TO_TEMPCOLV, imgname,first_caieq) + GLV_10_TEMPCO, imgname, first caleq); 
end; (delta_temp) 


(sunne 

procedure READ_tMAGE_IWFO( 
var 
ver 
ver 
ver 
var 
wer 
wer 


(procedure reads a Comm 


az 
el 


$f (oglirng © 2.0) then 
lope s# 0.009533 
else if (esl rng = $.0) then 
ep\_slope re 0.0207 
else if Cegl_rng # 10.0) then 
Slope 3= 0.04603 
if (egl rng = 20.0) then 
eq\_slope := 0.0871; 


GLV_TO_TEMP te 273.15 ¢ 
Cey|_slopetvaluee(ugi_temp-eg|_slopetegl giv) 3; 


writeln(output, ‘Ne eslibration equation available, check image name.‘); 


if unit_type © 4180! then 
GLY to_TEMP 1= isotherms * slope + intercept 


Temp 18 b / (in(¢a / isotherms) + ¢)); (Deg K> 
Temp te Temp - 273.15; 
Temp :9 Yemp * elope + intercept; (deg C) 
QLV_to_Tewp 12 Temp ¢ 273.158 (Deg kK) 


(Oeg Cd 


asnigntinfiie, ‘egl_cal.782’); 


while (mapensme <> eg{ name ) do 
readiniintile, egt_name, egl time, egltem, egiiglv 
es a 


(CRUDODUSUE TRADE U ERE POTTER SERES UCCs SASSER NURUEe REBAR eENEOEEEL -AnBE SNOT U REESE) 
+ Integer; 

1 image_name; 

first caleq : caleq_ ptr): real; 


AES UO CETNS CSUN SSRSE RSA RRARSeERS OEPEECU RE SHeRETEEE) 


image_name; 
Inveger; 
Integer; 
integer; 
integer; 
iriteger; 
real; 
real); 


delimited image information file, 
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AGAUTILS & OF 10 


Al5 


change to read apace del laited format after the OBU/VAL) AGAUTILS 9 OF 10 


ver 
infile t text; 
code fateger; 
ch,ehi, 

che : etring(t); 

Junko ecring(9); 
Junki0 + etring(10); 
Junk ateing(4); 


begin 
egslgn(inftle,nemer’. img! )t 
resetcinfiley; 
th smb; 
white (cher! ,/) do reed¢infile, en); 
eoad(intile, junké); 
vel ()Junké, year, code); 
read(intile,ch,cht,ch2); 
{f chet, then velceht , day, code) 
alae vel(chiech2 day, code); 
\f eh2et,’ then reedtinfile, ent) 
else read(infliie,ch, cht); 
valent month, cede)? 
reed(intlle,ch,oht,ch2); 
(¢ ch2e¢1? then waltcht hour, coda) 
else 
begin 
valchiech2, hour, code); 
read(infile,ch); H 
end; 
read( infile,cht,ch2); 
valcenterh2 minute, code); 
read( inf ile,ch,eh, junk?); 
val( Junk, a2 code); 
readcintite,eh, junki0); 
val( jumk10,e\, code); 
clesetinfile); 
end} 
(ansvevennaee 
Procedure READ _BASELINE_INFOC name | image nome; 
wer day : integer; 
var month : integer; 
yar year +: integer; 
var hour =: integer; 
var airate: integer; 
var et 3 real; 
ver el t peal); 
(oeseline image Information for DEN/VAL 13 Sept 90 bane! ining) 


Tritt TEE ELT LLL Ly 


ver 
err =o t integer; 
view +: Integer; 
begin 
day oe 13; 
month se 9; 
year 121990; 


vel (copy(name,3,2), hour err); 
minutere = 0; 
val (copy(name 5,2), view, ere); 
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el ga 


91,333; 


case view of 


1s 
2: 
33 
a 
St 
6 
T: 
6 
ay) 
103 
its 
42s 


ezr0157.0; 
0220159.9; 
0232 162.0; 
0210 166,5; 
0222167.0; 
02:9169,5; 
0z18172.0; 
0278190.0; 
0232192.5; 
az10195.0; 
azrn197,5; 
02:8200.0; 
arra202.5; 
begin 


0219235333; 
elie 92.333; 


begin 


(ARMAS OAnEee dues soeeenseuNeeCuerseESEReSER EUS UOURSUEAUrseusaeaenennannen) 


end. 
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WALT; 


AGAUTILS 


10 OF 10 
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A18 


Unit AgARange; 


interface 


wes 
printer, 
agaut{(s; 


cones 

hpin 782 » 0.456; 

vpin_782 © 0.456; 

hpin 670 © 0.3107 

vola_@70 © 0.310; 

rnerow = © 160; 

rrycol = ® $7) 

namin © 153 

target.width = 6; (eet width of target bon in meters) 


type 
ronge_image © array(0..engrow, 1..rnecol] of integer; 
uindowrec = recerd 
renege t integer; 
halae : integer; 
velre : integer; 
ord; 
windowarrey © arcay[0..numint of window_rec; 


ver 
renees t range 


(orevmcnnenersvonersenss Poerrres 


procedure read_range(dataloc +: string!2; 
imagerame : image_name); 


procedure find_bouncs_ype( = ax + real; 
al : coal; 
war lowbound +: integer; 
var upbound =: integer; 
{magename : image name); 


procedure find_bounds(ver lowound 1 integer; 
ver upbound 1 integer; 


imagename : imege_neme); 


function get_renge( row : integer; 
eol | integer; 
{magerione | image_name) 1 integer; 


procedure init_windows(ver windows 1: window arrey; 
\mageneme : imege_name)! 


procedure get_window(range integer; 
windows Window_erray; 

vee lowrange 1 Integer; 

ver uprange Inteyer; 
vor window window _ree); 


(oe tewancasremouncevaracessooenen 


imptementation 
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rari ey 
procedure read_range(dataloc : string!2; 
imagename : image_name); 


AGARANGE 2 OF 


i, | & Integer; 

infile ¢ file of integer; 
temp =o integer; 

apacnt : integer; 
etring2s + atring (23); 


begin 
(f copy(imagename,1,2) = ‘EG! then 
atring?S 16 dataloe © ‘—GL! + copy(imegename,7,2) ¢ '.rng! 
else 
{f copy imagename, 1,2) © ‘YG’ then 
atring2S s= dataloc + ‘YPG! + copy(imagename,7,2) + '. rng! 
else 
{f copy( imagename,1,2) = ‘AH! then 
string25 = dataloc + ‘APH! + copy((magenane,7,2) + ‘.eng! 
else 
{f copy¢imagename,1,2) © ‘AG’ then 
etring2?S 1» dataloc ¢ ‘APG! + copy( (mageneme,7,2) + ‘.rng’ 
else 
{f € copy(imagename, 1,2) * 'FD') and ( copyCimagename,11,1) = ‘'1/ ) then 
string?S :# dataloc ¢ ‘FID! + copy(imagename,7,2) 4 '1/ + ' Mme! 


else 
{f € copyCimagename, 1,2) © *FO') and ( copy. imagename, 11,1) = ‘2! ) then 
atring25 1@ detalec + ‘FID’ + copy(imagenmme,7,2) + 2’ 4 ‘.rng’ 
else 
{f copy( imagename,1,2) = HL! then 
atring2S := dateloc + ‘HTL! ¢ copy(imagename,7,2) 4 ‘.rng’; 


asulencinfite, etring2s); 
reset¢intile); 


if ¢ copy(imagename,10,1) © '7' ) then 
begin 
¢ For AGA7E2 ) 
for 1 1» 537 downto 3} do 
begin 
egacnt is trunc¢{ / (537/140) ); 
for } 18 1 to rrgcol do 
tead(infile, ranges(agacnt, j}) 


end 
end 
eles 
begin 
{ for AGA870 > 
for { 18 1 to 84 do 
for j te 1 to rrgcol do 
rend inf ile, rangest{,j)); 
for 1 :» 367 downto 1 00 
begin 
mgacnt 1 trume(l / (367/140) ); 
for j := 1 to rngcol do 
reed(intile, renges(agacnt, j)) 


end 
end; 
closecinfile) 
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end; 
AGARANGE 3 OF 6 
(ARRON SQRUOD SREY SAAREOR SSCS NERER CERN SERSERSNe REE SSAA E EEE Nee EE eed eKseEAEe) 
funetion get_renge(row 1 Integer; 
sol 1 Integer; 
imogername + image name) : integer; 


var 
newcol 1 integer; 
tmp_var : integer; 


begin 
{f ¢ copy( imegename, 10,1) © '7* ) then 
newol te trunc( 2+ (53 * col) / 140) ¢ For AGA7AZ > 
elee 
fewest t© trunc( 10 ¢ (37 " col) / 140 9; ¢ For AGAR7Q > 


if (reweel < 1) then 
feucol ta 1 

alee if (neweel » 57) then 
newcel so SP; 


get_range te ranges(reu, newcol); 
tep_ver 1" ranges(rew, newcel); 


(V ¢ © Coopy(lamgename, 1,2) © /¥G') or Ccopy( imagename,1,2) © “ygt) > 
end (tap var < 400) ) then 
Pet range:22000; (<--°FIX TO FILTER QUT GAD RANGE DATA FOR YP) 
end; 


CeneasveseReRessanasesyEeeeReansacensnensacteguens RUSAURGeaeeSeeeRneEEEEE) 


Precedure find _bounds ypg( a2 b real? 


el r real; 
ver \owbound 3 integer; 
var upoound +: integer; 
{magename : imege_ nome); 


(detorminas range limita of processed imagery to be applied 
for YUMA DEM/Val)> 

conat 
lower » 800; 
vpper » 29000; 


type 
date_pair © record 
ar: integer; 
el: real; 
end; 


ver 
(,)  t Integer; 
upper_done : boolesn; 
ridge t erray('..31) of data pair; 
bottom 3 Peal; 


begin 
ridgelt).ar ce 158 cidgett).el s2 92.5; 
ridgel2).az 15 140 ridgel2).el t= 92,4; 
chdye(3).an re 162 ridgell) .el s= 92.7; 
cidge(4) en 12 16 ridget4) el 1© 92.6; 
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tidge(S).az 1© 166 
Piden(6).an 1¢ 168 
tidget?] az 1s 170 
ridgel@].ax ta 172 
ridgel9) az te 174 


tiéee(10) aa 
eldget 1) aa 
ridgeti2) .az 
ridgelt3) az 
ridge (14) .a2 
ridgeltS) az 
ridgelté) az 
tidgel17) .a2 
ridge (18) aa 
ridge 19) aa 
Fidge (20) az 
ridge (21) a2 
ridget22) ,az 
tidge (23) .a2 
Cidew (24) az 
idge (25) az 
ridge (26) a2 
Pidge (27) .az 
ridge (26) .a2 
ridge (2?) a8 
Pidgel30) .az 
ridge (31) .a2 


upbound 1° 1; 
upper tiona is falee; 


Vie ty 


sate 
10178 
18180 
19182 
10th 
18t8é 
10188 
te 90 
inf9a 
10196 
10406 
19198 
18200 
tez0a 
10204 
10206 
1aade 
10210 
102i2 
taza 
te2he 
102346 


ridge(S) .e( 

ridge(6).e4 

ridge (7) .al 

ridge (8).e\ 

ridget9) .al 

ridgeCi0].¢1 
cidgeti1).et 
eldge(12).ei 
ridge(!3).¢4 
ridge l14).ei 
ridge (15).¢1 
ridge (16) .el 
eidee(17).¢1 
ridge (18) .01 
oidge (19) el 
ridge (20) .01 
ridge (21) .el 
eidge(22).e\ 
edge (a3) .e1 
eidge (24) el 
ridge(23) .e1 
lege (26) 01 
ridget27) 01 
ridge {28} .el 
fridge (2) .@| 
idpe (30) .@1 
eidgel31) el 


OS 


18 92.6) 
88 92.6; 
as 92.5; 
15 92.5) 
12 92.5; 
0092-4) 
rege. ty 
1091.8; 
1091.8; 
oO ey 
reV1.7) 
1091.8) 
1091.9) 
1092.9) 
1992.9) 
1091.9; 
109204; 
2992.0; 
1991.0; 
1091.9) 
1091.7) 
Ar 
1910; 
1091.0; 
191,93 
1997.0; 
1v97.0; 


Wille ridgeci) agcround(er) do f(ieiet; 
battem re ridgeti).et; 
Lire dp ote PO; 
while net(upper_dene) do 


begin 


(f GET RANGECY , j, imagename) «< Upper then 
begin 


upbound 1° {; 
upper_done t= true; 


end) 


Pte; 


ond; 


if ( copy(imagename, 10,1) © '7' ) then 


Vowbound 16 round( (bot tem-@19*¢1,0/¢(vpix_782/1000,09"57.3))) + 70 


\uwbound 18 round( (bat tom-et )°¢1,0/((vpix 870/1000.0)"57.9))) + 70; 


(f Upbound<10 then upbound :@ 10, 
{f lewbound>120 then luwbound 12 120; 


end} 


(Ramana Petes naeNe Resa ans DURE H ea LSU AAS Ee RESeAQueseyRAeeENOcRRASAMEeReReE AE) 


procedure find_boundsa(ver loubound 


comet 
lower = 300; 
upper = 29000; 


ver Upbound 


1 integer; 
1 integer; 
\magename : image_name); 
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A2i1 


ver AGARANGE 5 OF 6 
is] : integer; : 
done 1 boolean; 


begin 
dene tm PALEE; 
bre tg y 18 PO; 
walle netiaene) do 
begin 
(¢ cer nanan! , J, lmagename) < wpper then 
begin 
Upbeund 18 1) 
done 1s true) 
ond; 
fries} 
ord; 


dene 16 PALGE; 
t 10 140; J 18 70) 
while motidene) de 
bepin 
if Gat mavamd! , J, \mageneme) > lower then 
begin 
\evbeund 16 1) 
dene 18 true; 
ond; 
liete ty 
ond; 


1 Upbound « 10 then upboied r= 10) 
1f lewbound » 120 then \ewbound 1# 120; 
ond) 
(HRCRDP DORN ERAAR PRADA DEARORSACORARARASAEGeSSRaAeR RAs AaeE A BaLReeeeennenAneD) 
presedure Init windows(ver windows 1 Window_erray: 
imagenema ¢ image name); 


ver 
4 + integer; 
begin 
for | 1 mamin dounte 1 do 
with windows (numwin > 1) do 
begin 
heise re 2 Oe io 1 
If CCtruned( 4 2) * 2) © trunedi / 2)) then 
walre ie lof 
else 
woiee im fy 
{ft ( eopyCinagenam, 10,1) © 'f’ ) then 
Forge 12 trunc(terget_width / einthsizse * (hpin_782 / 1000))) 


else 
renge ie trunc(target width / ain¢hstze * (hpla_@70 / 1600))); 
end) 
ed; 
(passceagces: OTST CDARTOSUANGRETEESS DUN ESSRSERUSOSHRSUSECSDEAAHeDELeBEeeS) 

procedure get_window( range t integer; 
windows : window array: 

var lowrange « integer; 

ver Uprenge 1 integer; 
war window 1 window rec); 
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A24 


Unit aAGAMetre; 
interface 


ees 
aseutile, 
eparange, 
epetgtin; 


eorat 
histostze = 256; 


type 
hisyopte = “histogram; 
histogram = record 
bine t array(-235,.456) of 
end; 


integer; 


{ate emacs eemeneesensnensaseaneneereceasenesseusecmecsnnonscgus) 


Procedure sentraat(var hiete 
leweound 
*upibouind 
\mage 
wlndews 
ver cent_img 
firet_tgt 
{mageneme 


Procedure sontrast2(var hiete 
\owbound 
‘upbound 
image 
windows 


var cont_iag 


t hieteptr; 

+ inteser; 

t Integer; 

1 (mage_array; 
1 window array; 
t (eage_array; 
t target ptr; 
t teageneme); 


t hleteptr; 

+ integer) 

t integer) 

t image_array; 
t Window array; 
1 cont_array; 


{magename 1 image name); 


function GTPCeontrest 1 real; 


hista 1 histoptr): real; 


function GLY_GTP(hot t Integers 
hist 1 QLV histo) 1 real; 


function PERCENTILE( Ver histo 
percent 


[ Coe ee ee 


tap ementet ion 


1 hietoptr; 
1 realy: integer; 


(Cenennecenneopaner erste sene Ue sULet canes nan eeaees Basen emeneedeesuaEs) 


precedure contrast(ver hiete 
(owbound 
upbound 
image 
wlindows 
var cont, lag 
first tot 


type 
tet window © record 


histoptr; 
Integer; 
integer; 
\mage_array; 
window _erray; 
image_array; 
tatget_ptr; 
image_name); 
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integer; 
integer; 


xs 
y} 
x2: Intege 
y2: integer; 
and; 
tgtlarray * eremy(!..20) of tot_window; 
tot_itee # “tet, (let_rec; 
tat_iist_rec © record 
mum + integer; 
mont 1 tet list; 
end; 
(eres eenmeene PPPPITTT TTT reer rrr 
procedure make_tgt_windows( firet tgt 3 targetptr; 
var target_windowa : tgt array} 
var num_tete t integer); 


ver 
curr_ptr + target_ptr; 
i 2 integer; 


begin 
cure ptr ce first tat; 
ite 43 
while (curr ptr*.next_tgt @ nil) do 


vesin 
target windows (1].x1 18 curr_ptr*.target.centerx - trunc(eurr ptr*.tgtuind.hsize / 2); 


target_windows(i).x2 t@ curr_ptr*.target.centern ¢ trunc(curr ptr’. tgtwind. helze / 2); 
target_windowsil].y1 t* curr_ptr®. target.centery - trune(curr_ptr*, tatwind.veize / 2); 
target_windows{i].y2 t= curr_ptr*. terget.centery ¢ trunctcurr_ptr*. tytwind.veize / 2); 
curroptr 18 curc_ptr*.neat tot; 
brelet? 

end; 

ram_tgta re d+ 1; 

ond; 
PPryererrrerrerirrrrr irre rir re rr iisS | 


procedure sort_tergeta(ver terget windows : tgt array; 


rum_tgts 1 integer); 
var 
temp_terg t tgt_window; 
i, J : integer; 
begin 


for | t© (ram tats © 1) dowta 1 do 
for j r= 1 to 1 da 
{f (earget_windows(/].y2 > terget_windows[) + 1).y2) then 
begin 
tump_targ 1 target, windowslj ¢ 12; 
target_windows() + 1) 18 target_windows{)]; 
target_windows(J} :¢ temp_serg; 


procedure check target(row 1 Integer; 
window : window rec; 
vor tgtliat : tet_list: 
ram tte 1 integer; 
target windows : tet_arrey); 
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{ 1 integer; 

top 2 integer; 

bottom 3: integer; 

currtgt : tetli let; 
begin 

if Ctetlist <> nil) then 

begin 


currtgt :* tytlist; 

while (eurrtgt «> nit) do 

begin 
tetlist s© currtgt*.next; 
Cispose(currtgt); 
currtgt ss tetlist; 

ord; 

end; 

new(totl iat); 

currtgt :@ tetlist; 
currtgt*.next re nil; 

{a} 

tap r+ row + window. veize; 
bottom :e row + window. veize; 
tor { 18 1 to mum_tets do 
begin 

if (Crarget_windows(}].y1 >= top) and (terget_windows (I) .y1 <* bottom) ther 

begin 
neu(currtet* next); 
currtgt 18 currtgt*. next; 
currtgt*.num te i; 
curetet*.mext re nil; 

end; 

{f (Ctarget_windows (i].y2 >© top) and (terget_windowsli}.y2 <e boxtom)) then 
if (eurrtgt*.mum <> i) then 
begin 

meutcurrtgt* ment); 
euretgt t@ currtgt*.next; 
currtgt* rum te 4; 
currtgt*. next re nil; 


procedure update targets(var target windows : tgt_array; 


ver mum tots : integer; 
window > window rec; 
tow : integer); 


ver 
1, J t integer; 
top : integer; 
begin 
top 5 tow > window.vsize; 
for | re 1 to ram tgts do 
begin 
if (terget_windowsii).7) ° top) then 
begin 
tor j :8 1 tou mum_tgts © 1 do 
terget_windows(j) :* target_windows {J 
mam _tgta (© ram tgit «1; 
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+ integer; 
rum_tots + integer; 
tgtliet + tet lise; 
target_windows : tgt_array; 
accuminner t longint; 
accumouter t longint; 
rampixel t integer; 
@eaninner t real; 
meancuter + teal; 
AGA_Contrast : integer; 
tenge 3 integer; 
window : window_rec; 
lourange : integer; 

. uprange : integer; 
currtgt t tetllist; 
procede + boolean; 
top + integer; 
bottom t integer; 
lefe t integer; 
right t integer; 
current_tyt + integer; 


begin 
for { :# +255 to 256 do 
histo*.bina(i; := 0; 
make_tgt_uindows(firet _tgt, varget_windows, mum_tgte); 
sort_targete(target_windows, tmm_tgts); 
range t# get range(upbound, left_junk, imagensme); 
get_Window(range, Windows, lowrenge, uprange, Window); 
, k t= Upbound + window. velze; 
Sethist s= nil; 
for k += upbound to lowbound do 
begin 
check tary’ t(k, window, tgtlist, rum tots, target_windows); 
range := get range(k, (left junk * window.hsize © 1), tmagename); 
if (range < (orange) or (range > uprange) then 
pet_window{range, windows, lowrange, uprange, window); 
for 1 t@ (left junk + wWindow.nsize + 1) to (mumcol > Fight junk + window.nsize) do 
begin 
procede := true; 
currtet := tgtlist; 
left :¢ L + window. hsize; 
tight :2 i ¢ window.halse; 
top 1 k + window. veize; 
bottom := k + window.veize; 
while (curctgt*.mext <> nil) do (identify positions where contrast box overlaps terpets) 
begin 
current tpt 3 currtgt® num; 
if (Ctarget_windowslcurrent_tgt).x2 >© left) and (terget_windowsicurrent_tot).a2 <# right)? then 
procede := false; 
it (Ctarget_windowslcurrent_tgt).x1 >= left) and (terget_windowslcurrent_tat).al se right)) then 
procede :* faise; 
if (Cterget_windows [current _tgt).yl 22 top) and (carget_windowslcurrent_tgt].y! <= bottom)) then 
procede := faise; 
{ft (Cterget windowslcurrent tgt).y2 >= 
procede 1° fatse; 
currtgt := currtgt*.mext 


top) and (target _windowslcurrent_tgt).y2 «* bottom)) then 
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end; AGAMETRIC § OF 8 
if procede then 


begin 
accumouter :@ 0; 
rampinel se 0; 
renge 1° géet_rengetk, L, imegename); 
{f (range < lewrenge) er (range > uprange) then 
#et_window(range, windows, lowrenge, uprange, window); 
left t2 (= wlndow Asie; 
right te 1 + window. hsiae; 
top 12 k * window. vsize: 
bottom 18 k + window. veize; 
for | te Ch - window.vaize) to (k + window.vsize) do 
for ] t= left to right do 
begin 
writeln(output, (he, 1:5, ¢ he 4, ke, 6 fot, ped, Lee, US, § lete m , bettss, 
‘ eight # ', rights5); 


accumouter t= sccumouter + image{i)*.datalj); 
tumplxel t= mumpixel + 1; 
end; 
accusinner re 0; 
for { te Ck > trune(window.wsize / 2)) te (k © trunc(window.vsize / 2)) do 
for J] r= (i > trunc(window.neize / 2)) to (1 + trunc(window.msize / 2)) do 
accuminner 1@ accuminner + immgeti)*.datel)); 
eccumouter i# accumouter - eceuminner; 
rumpixel t+ numpixel + (window.velse * windowisize)s 
manninner t* accuminner / (window, veize * window.hsite); 
@eanouter 3" accumouter / mumpixel; 
AGA contrast 1s round(meanouter » meaninner): 
histo*. bins (AGA_cantrast) 1© histe*.bins (AGA_contrest) + 1; 
hinto*.binel2s6) re histo*.binst2S4) + 1, 
ord; 
end; 
Update targets target_windows, mumtyts, window, k); 


(tec etegusesacenscununeassess MSGRUERnevesaeecceeeeseEsuEEES) 
Procedure contrest2(var histo : Wistoptr; 
Vowbound integer; 
upbound integer; 
image : image array; 
witdows window_erray; 
cont_img : cont_arrey; 
imagename : lecce_name); 


var 
thd, uy a integer; 
eccuminner longint; 
eccumouter longint; 
rampine| : integer; 
@ueninner real; 
m@eanouter reel; 
AZA_Contrast integer; 
ronge integer; 
tab recht : window rec; 
lowrange t integer; 
uprenge : integer; 
Procede + boolean: 
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t t integer; 
eric 1 integer; AGAMETRIC 6 OF 8 
left s integer; 
right t integer; 
irmer_pixela +s integer; 


begin 
for i ss 1 te numrow do 
Rew( cont lag (i); 
for { ta 1 te mumrow do 
for k :¥ 1 to mumcol do 
cont_imet{)*.data(k) 2* -$28; 
for 1 :® «253 to 284 do 
histo*.bineti) 28 0; 
range 1 get_rangecupbound, left_junk, {magename); 
Ott window(range, windows, lowrange, Uprenge, window); 
k 18 Upbound + window, valze; 
for k 1® upbound to lowbound do 
bagin 
Fange 18 get_rangeck, (left_junk + window hslze + 1), imagensme); 
{f (range < (owratige) or (range > uprange) then 
Wtt_window(range, windows, lowrenge, uprange, window); 
for { 10 Cleft junk © window.helae + 1) to (numcol +» right junk - window. haize) do 
begin 
left 101 + window hatze; 
eight 12 | © window.hedze; 
top :# k + window. velse; 
bottom :0 k ¢ window. valze; 
brocede t« true; (process all pixels) 
{f procede then 
begin 
accumouter 35 0; 
mnpixel i» 0; 
Inner pixele :# 0; 
range ty get_range(k, |, lmagename); 
{f (range « lowrange) or (range > Vprange) then 
Ott window range, windows, lowrange, upranga, window); 
(eft i= + window. bsize; 
right te 1 + window. hvize; 
top 18 & + window. vsize; 
bottom := k + window. veize; 
for f rs Ck + window, vBize) to (k + window.wsize) do 
for | :© left to right do 
begin 
writelnfoutput, ‘6a ', 1:5, Ske ', hed, ¢ 224, Jb, bt, lib, 4 bette ee, left:5, 
* right =", right:5); 


accumouter 1s accumouter ¢ imageli)*.date(j); 
cumpinel :2 numpinel + 1; 
end; 
accuminner :* 0; 
for t rm (Ck + trunc(window.veize / 2)) to (k © trunc(window.velze / 2)) do 
for } :e Cl - trunc(wirndow heize / 2)) to Cl + trunctwindow.heize / 2)) do 
begin 
inner pixels i© inner_pinels ¢ 1; 
accuminner := accuminner + image(i)*.detel)); 
ond; 
eccumouter :* accumouter - eccuminner; 
numpinel :* mampixel + inner_pinels: 
@eaninner := accuainner / inner_pixels; 
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meanouter := eccumouter / rumpixel; 

ABA_contrast 25 round(meaninner: meancuter); 

cont_img (kJ*.data{l} :© AGA contrast; 
hiete*. bins tAGa_eentrast) :© histe*.bins{AGA_contrast) + 1; 
hiato*.bine (236) te histe*.bina (286) * 1; 


Bernas Tritt iit hy 
function GYP(contrast : real; 
hieto t htstopte) 1 real; 


var 
eccun =: longint; 
teap t real; 
{ t integer; 


begin 
eaccum t= 0; 
far | s= +235 te (round(eentrest)-1) do 
eccum 1@ eceum + histo*.binsti); 
temp :© eccum / histo’. bins (256); 
GTP ie temp; 
end; 


tunetion GLV_GTPchot t integer; 
hiet 1 GLV_Aiete) : reel; 
var 
accum =: longint; 
temp =: real; 
i : {nteger; 


vein 
eccum ru 0; 
for i t= 0 to hot do 
eccum t= accu > hiet({); 
temo :© eccum / hint (256); 
GLV_OTP i= temp; 
end; 


(cgenesceseenarenuses reset cess BEGeDESHSSRSCURSASAUSESENS SGRS RA PERORERESEEEEEEEEEeRERE) 


COLT EE CL Oo ey 


A30 


fumetion PERCENTILE(var histo : histoptr; 
percent : real): integer; 
ver 
stop_at : inteser; 
‘ 1 integer; 
eceum =| longint; 
beain 
ttep_et 1© reund(percent* int¢hinto*.bins{256) )); 
eccum = 18 0} 
1 is -835; 
while Caccum<stap_at) do 
begin 
eccum := eccum ¢ histo’. binali); 
{ried 
and; 
PERCENTILE se 1; 
end; (end percentile) 


Appendix A 


Metrics image Processing Software Source Code 


(enconmqcevavsssuesanasecevannatnunensesesnanexsneneeuEsneaneucaucensennescucassuneen) AGAMETRIC 8 OF 8 


end. 
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ne a ct ERO rin + 


unit egaceapt; 
interface 


uses 
, 
ere 
agarange, 
eqnut ( (5; 


conet 
angle_782 = 0.000456; 
angie,870 = 0,000310; 
Gia = 14.0; (Colze of clutter box in meters) 


type 
byte_earcier © array(0..10000) of integer; 
real_carrler @ array(0..1000) of real; 


function VARIANCE(dete : byte carriers! real; 


precedure senp_reynid( {ange 1 {mage ervey; 
(megerame + image name; 
upbound =: integer; 
leowmeund ¢ integer; 
ver clutter 3 real; 
var reynold : real; 
std_dev +: reat); 
f CO errs ween Cee | 
lnplementat | en 
(ebncens Ocenas denen cdueseus eben SER ESR EREEEBURSSBAeaaseENecEsESRaTERee) 
function VARIANCE(date ! byte_cerrier): real; 


+ integer; 
1 real; 
> reals 
t real; 


asm 0.0; 
for | 18 1 to datal0) do eo 18 # + detel}); 
ove te a / dace(Ol; 
aver :2 0.0; 
for j 12 7 to data(0) do 
begin 
p t® sqriratalj)-ave); 
Sver 12 Svar + p; 
ore; 
VARIANCE se swar / (detal0)+1); 
ord; ( end VARIANCE > 


(eenes seen es eneemes Coen erence cane anne ene eeweenenssneceneaann ateeeeners } 


procedure comp_reynid¢ image (mage_array; 
{mogename : (mage name; 
upbound integer; 
lowbound integer; 
var clutter teal; 
ver reynold : real; 
Btd_dev real); 
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var 
byte_dat 1 byte_carrier; 


AGACOMPT 2 OF 3 


ree! _dat 1 reel carrier; 
slre_pix 1 real; 
wldeh 1 integer; 
vert_boxes 1 integer; 
horiz_boxes §=«1 integer; 
residual 1 Integer; 
row start + integer; 
col _etart t integer; 
iy), fb, Jb,k kbs integer; 

t reat; 

t (nteger; 


(based on center image range? 
Prange 29 get_renge(trune(( lowbouncdtupbound)//'2), trunc((numcoleright Junk)/2), imagename); 


1€ ¢ copy(imagename, 10, 1) = ‘7! ) then 
size_pixn ie range * ein angie_782 > (nominal pixel site in meters) 
else 
alzepix te range * sing angle_B870 >; (nominal pixel size in meters) 
uldth em trunc(dia / elzepix); (pixel dimension of unit box) 
vert_boxes © (lowbound-upbound) div width; (number ef bones in vertical directions) 
Cen_stert ¢= Upbourd; 
horiz bones :© (nuncoleright_Junk) div width; (rumber of boxes in vertical direction) 
col_start :s t; 
1f (vert_boxes<>0) end (horiz_boxes<r0) THEN 
begin 
Ccompute and store variance of pach box) 
k re 0; {box counter) 
for i = 1 to vert boxes do (boxey in vertical) 
for j :@ 1 to horiz_boxes do (ooxes Im tortzontat) 
beatin 
bb 18 0) (pixel counter in box) 
Cleed carrier array and compute box variance) 
for ib :© (¢i+)*widtherow atart) to Ci*width+row wtert+t) do 
for jb 1m (Cj-t)*widthecol start) to (J*width>1) do 
begin 
kb re hb; 
byte_dat (hb) r= imagelib)*.datat jb); 
end} 
ktek et: 
byte det{O) :© hb; (pinele in box) 
real atk) :* VARIANCECbyte_det); 
end; 


accua 18 0; 
real_dat(O) ce int¢k); 
for k t# 1 to round(reai_det(0]) do accum 18 occum + real datth); 
clutter re eqrt(eccunvrent_dat(0)); (convert to temperature) 
If atd_dev «> 0 then 
reynold r= (atd dev - clutter)/ std dev 
else 
reynold 15 6858.5 
end 
vise 
begin 
clutter := 9999.9; 
reynold :¢ 9999.9 
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end; AGACOMPT 3 OF 3 


(onscensouccenguconsunsasenauacreneccauazadenernsennsecnmnsSnaeneurcen) 


ond, 
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PROGRAM: VISCCD.PAS 
MAIN PROGRAM FOR 
VISIBLE IMAGES 


READ FIRST 
SITE [O BE 
PROCESSED 


BE PROCESSED 


READ NEXT 
SITE TO BE 


CODRANGE,. PAS 


CCDUTILS.PAS 
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A36 


COMPUTE MORE 
IMAGE 
METRICS 


CONVERT IMAGE 
METRICS FROM 
GRAY LEVEL VALUES 
© TEMPER 


TURE 


CCDUTILS.PAS 


CCOCOMPT.PAS 


CCOMETA.PAS 
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unit CCOUTILS; 


interface 


uses 
dos, 
ort, 
tpextnem; 


conat 


type 


ond; 
image_array 
work, ree 
image_row 
videolartay 
video, row 


end; 


end; 


C ROTES This program normalize GLv to exposure timesd.! (If exposure time ie 
different from 0.1)> 


huarow wy; 
murcol = 576; 
feft_junk = 3; 
right junk # $7 


Qlvptr = * “*OLV; 
Gy o= record 
bine 1 arvay(0..16300) of lengint: 


© array(\,.mamrew) of hugeptr; 
© arcayi!,.1024) of integer; 
= arreyl!,.mumcol) of integer; 
© array{\,.mmrow) of hugeptr; 
@ array[t,.102¢] of byte; 


comment info = erray{!..2) of char) 
image_comment © record 

zoom =} comment_ into; 

fatop : comment info; 

filter + omment_ into; 

June ss errayli..94) of cher; 


image_pos = record 
col t integer; 
row : integer; 


fmage_header * record 


comment ¢ (mage comment; 
eaptimel 1 word; 
exptime2 1: word; 
« exp_time 1 longint:) 

Juna! 1 atrey{l,.6) of byte; 
img_time ot longint; 
jeg_org ot imege_pos; 
img teng : tmage_pos; 
fap_bin 1 (mage_pos; 
gain 1 integer; 
read_rate : integer; 
Junk? 1 atray(),.16) of byte; 
num_frome 1 integer; 
jog type 1 integer; 
jumk$ 1 array(!,.8] of byte; 

end; 

image_name « string (8); 

atring!2 «© atring(12); 

atringié = = atring (16); 

q a rrr rrr rr eree rt rer ee eee ee >? 
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A37 


ver af 
wer el 


Procedure read_image(var image 


ver Givhist 
var min 
ver man 
var percds 
ver perc0S 
ver eadien 
var mode 
ver meen 
ver g00m 
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Procedure racd_image(var image t imege_array; CCDUTILS 2 OF 8 
nme t image_name; 
var GLY¥niet +: GLYptr; 
var ain t integer; 
ver max : integer: 
var perc#S =: integer; 
var pereoS : integer; 
var wedian =| integer; 
var mode t integer; 
ver ween t real; 
ver soon 1 comment info; 
ver fetop t comment_info; 
ver exptim ¢ real; 
var filter 2 comment_info; 
var etart 2 hugepte; 
var entropy : real; 
loubound : integer; 
upbound : Integer; 
Gata_loc 2 stringté); 
Precedure create di splay( image 1 image_array; 
ver display 1 video array; 
ver etd_dev + real; 
ain t Integer; 
ax t integer; 
moon t real; 
var skewness : reel; 
ver kurtosis 1: real; 
lowbound : integar; 
upbound : integer); 
Procedure READ IMAGE IWFO( name 1 image name; 
ver day «1 inteyer; 
ver month «integer; 
var year 1 integer; 
var hour 1 integer; 
ver minute: integer; 
var at 1 real; 
var el t real); 
Procedure READ_GASELINE INFO( name : imege_name; 
war day: integer; 
var month 3 integer; 
var year =: integer; 
war hour: Integer; 


var minute: integer; 


Co nenartneaersccrceessrees Peerrrrrert rrr erert rere rrer rr rrerirrs) 


OUORNs MNEs eRe eOt COONS CE ESEBESET Ages eeEEseeeneasE) 


: real; 
: teal); 


image_arrey; 
image_name; 
GL¥pte; 
integer; 
integer; 
integer; 
integer; 
integer; 
integer; 
real; 
comment_into; 
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wer fatop : 

ver exptia =: 

ver filter 1 comment_info; 

var atert + hugeptr; 

var entropy : real; 
: 
: 
t 


comment into; CCHUTILS 3 OF 8 


real; 


lowbound : integer; 
upbound integer; 
data_toc 2 string!6); 


(rneeeccenee eenenaen PrerT errr rrr rrr rr rrr rrrerr irr rT Serr triers | 
procedure read_hesder(date_loc : string!6; 
name =: image_name; 
ver toom = 3 comment_into; 
var filter t comment_info; 
var exptia t real; 
ver fatop 1: comment_info); 


vir 
infile ¢ file of imape_header; 
header : imege header; 
begin 
essigncintile, data_loc ¢ name + /.cot' D7 
reset(infile); 
resd(infile, header); 
cloae(infile); 
zoom :© header. comment .200m; 
filter 1 header, comment. filter; 
fatop ra header .comment. fstop; 
header. axptine2? t= suap(header,exptimez); 
exptim 1" header,exptime2 / 10; 
end; 
(ebsvnenavensecnmccesacccenecctousssusarwarrasinesccscmsscceser) 
function GLV_PRRCENTILECvar hiete ¢ GLYptr; 
ni real; 
percent : real)iinteger; 
var 
@topat 1 real; 
( + integer; 
ancun: real; 
begin 
etop_et 16 percent*n; 
accum = =o1# 0; 
{ 1m 0; 
while Ceccumestop_ax) do 
begin 
accum t= eccum + histo*.cinell); 
{re be t; 
ord; 
GLV_PERCENTILE :© 1; 


Beemer eesecnesemascaecene} 


intitle file of integer; 

atorerow errey[0..(rumcot - 1)) of integer; 
work : work row; 

gridsize =: ward; 

oridsizek : word; 

1, J 1 integer; 
accu 1 real; 
percentil! 1 teal; 
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aiddie real; CCDUTILS 4 OF & 
percentil2 1 real; : 
Junk integer; 

fageize word; 

done beolean; 

dane! boolean; 

dore2 boelean; 

preb real; 

percentil real; 

totelrum real; 

temp integer; 

tespptr hupgepter; 

eddrece lengint; 

tecter integer; 

freq max longint; 


begin 
Pead_hesder(dete_lec, mame, room, fliter, exptia, fetop); 
factor t= reund(expt ia/0. 10); Cexpcaure tise »ormal ization) 
sesioncintiie, date_loc + name + ‘.ced’); 
reaetcintile); 
for 1 89 0 fo 14300 do Givhist*.binsf{) 1= 0; 
far 1 se t te 80 do read(infile, Junk); 
eridsize © ¢izeot(wark); 
eridsinek t@ trunc(erideize / 1024); 
eridaize t= round(gridsize / 2); 
warkectaum start); 
for i 32 1 to mmrow do 
begin 
image (|) :© getextmemorideizek): 
{¢ Cimapeli] © mit) then writein(output, ‘OUT oF WAN’); 
end; 
eceum 18 0; 
Man 1® +2000; 
win te 20000; 
for j := mumcol downto 1 do 
for 1 18 1) te eumrew do 
begin 
reed(infile, teap); 
temp :* (teep div tector); (normalize GLV to exposure times0.1) 
4f (teap> 16209) then temp: 216299; 
1f (temp<0) then temp: e0; 
4f CC) > Lett junk) end ¢j < rumcol - right junk) end (1 <= lowbound) end (i >= upbound)) then 
begin 
f (temp > mex) then max :* temp, 
tf (temp < min) then min te temp; 
eccum :© eccum ¢ temp; 
GLvniuc".bineltemp) 10 GLVi let*.bineltemp) + Vy 
Gi Vniet* bins 116300) :© Gtvnist*.bing (14300) « VY 
ond; 
amveentnom iangeli}, ptrtohuse(addr(work)), eriduize); 
work(j) = temp: 
Govees tam ptrtohuge(eddr(work)), image(ti, arideize): 


address :© longint( iaage!')); 
tenpptr re hupeptr(eddress + (j - 1) * 2); 
moves tmem pir tehuge(edde(temp)), tampptr, 1); 
ond; 
close(infile); 
Maen :¢ sccum / Givhist*. bine 116300) ; 
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perc0S :© GLV_PERCENTILE(GLVhist, int(GLvnist*. bine (16300) ),0.05); 
wedian 16 GLV_PERCENTILECGLVh ist, int(CLYhist*.bing (16300) ),0.50); 
percds GLV_PERCENTILE(GLVh Ist, int(Glvhist*. bins (16300) ),0.95); 
entropy: 90; 
treqmax 19 0; 
for iremin to mex do 
begin 
If ¢ givhist*.bing(1] > freqmex ) then 
begin 
frecq_max 1= givhist*.binsli}; 
mode 15 | 
end; 
prob re GLVhist*.bins li) /int¢GLvhist*. bins (16300) >; 
1f (prob > 0) then entropy := entropy * (prob * In(prob)); 
and; 
entropy <* -entropy; 
ord; 
(em mereseneransenese suns es BURG CORE SenERee Ree Re tees Nene Ree ERE EE RR eRe se Meare aEKsseseee) 
procedure create_display< lmage image_array; 
var display video_array; 
ver etd_dev real; 
ain integer; 
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mean 

ver skewness 

war kurtosis 
lowbound : integer; 
‘upbeund integer); 

r¢ wer erew eee Preece eee eee eee? eowetervesrre Pere ee eee eee ee ) 
tunetion INT _TO_GL¥(value 1 Inte 
ain 1 integer; 
max 1 Integer) : byte; 


Peal; 
real; 
real; 


t 
t 
t 
t 
ann t Integer; 
: 
1 
: 
: 


var 
temp : real; 

begin 
tempt e((vaiue-min)/(mancmin)); 
1€ (tempt) then temps; 
1 (Ctemp<0) then temp: 50; 
int_tolgly = trunc(temp’63); 


: video tom; 
1, j : integer; 
eridsize =: word; 
videoaize ; word; 
eridsizek : word; 
eccunt real; 
eccum2 real; 
accumd : teal; 
equere 2 teal; 
tewp s real; 
rappixel ; real; 


begin 
uridstae := eizeof(work); 
Orideize s+ trumc(gridsize /2); 
videosize :© eiteof(video); 
Eridsluek :¢ trunc(videosize / 1026); 
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videosize s= trunc(videosiza / 2); CCDUTILS 6 OF 8 
eccual 32 0; 
accu? 32 0; 
accucd 10 0; 
rumpixel := 0; 
for 4 :# 1 to 1026 do videot{) se 0; 
for { 28 1 to numrow do 
begin 
diepley{l) 2% getextmemgridsizek); 
{f (display(() © nil) then writein(output, “GUT OF RAM (CREATE DISPLAYD’); 
moveextmen( image it), p' rtohugeceddr(work)), gridslze); 
for j 1 1 to mumeol do 
begin 
videot)) 12 int_to slv(work(j}, min, max); : 
if CJ) > deft junk) end (J < mumcol > right_junk) and (i c= lewbound) end (i »* upbound)) then 
begin 
rampixel s@ nunpixel + 1; 
temp re work{j) » mean; 
equere t= sqr( temp 
eccus! := eccunl ¢ square; 
eccum2 16 eccumt + (square * temp); 
eceuaS t= eccumd + (square * square); 
ond; 
end; 
moveentmanptrtohuge(adde(video)), displayl{), videosize); 
end; 
eccuat := accusl / (rumpixel - 1); 
Btd_dev re sqrt(aceunl); 
eccum2 12 accum2 | (accuml © std_dev); 
skewess :© accune / mumpixel; 
accum3 re sccum3 / (eqr(accum!)); 
kurtosis :# (accumd / mumpixel) - 3; 


end 


Ranecersecan sucunensensannuse) 


procedure READ IMAGE_INFOC name. image_name; 

var dey : integer; 
var month : integer; 
var year +: integer; 
var hour ; integer; 
var minute: integer; 
var at : real; 

var el ; real); 


(procedure reads @ comma delimited image infcrmation file, 
change to read space delimited format after the DEM/VAL) 


ver 

infile : tent; 
code =: integer; 
ch,ent, 

ch2 : astring({1; 
jurnkO os atring(9); 
Junki : etring(10); 
jureé st string(4); 


begin 

spsigntinfite,namee’, img’); 
Peset(infile); 

ch vel’be; 

while (cho’,') do read(infile,ch); 
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read(infile, Junks); 
val(junké, year ,code); 
read inflle,ch,chi,ch2); 
if enae',? then val(eht day, code) 
elae vel(chiech2,day,code); 
if chaz’, then read(infile,ch1) 
else read(infile,ch,ch4); 
vaitch! month, cade); 
reed(infile,eh,cht,ch2); 
if chéerr’ then valceh! ,hour,code) 
else 
begin 
val(chi+ch2, hour, code); 
read(infile,ch); 
end; 
read infile,chi,ch2); 
vai(chlech2,minute,code); 
read( infile,ch,ch, Junk?); 
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vai (junk¥,az, code); 
reed infile,ch, junk10); 
val(junki0,e(, code); 
closatinfile); 
end; 
(espe cannysvsnsaarseeensnannmnnensecunssesennsnaeetneensneneanesnyeecene) 
Procedure READ_BASELINE_INFO( name : image_neme; 
war day 1 integer; 
var month : integer; 
ver year : integer; 
var hour +: integer; 


var minute: integer; 
var az t reat; 
var el 2 real); 
(baseline image information for DEM/VAL 13 Sept 90 basel ining) 
var 
err + integer; 
view : integer; 


begin 

day <=» 13; 

month 12 9; 
year 191990; 
val (copy(name,3,2),hour,err); 
moinute;s 0; 
val(copy(name,5,2),view,ert); 
et :2 OL.333; 
case view of 

1; axse1$7.0; 

2: oxt0159.5; 

3: 0210162.0; 

&: on18166.5; 

5: 02:%167,0; 

6: a2:2149.5; 

7: ansei72.0; 

O: 02:0190.0; 

9: aus0192.5; 

VO: 02:9195.0; 

VWs ans897.5; 

V2: on:8200.0; 

13: a2:©202.5; 

14: begin 
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a21eZ3$ 333; 
elre 92.333; CCDUTILS 8 OF 8 
ond; 
elee 
begin 
writein(/incorrect view @ read from baseline file \iat'); 
WALT; 


(euaeeesonccanenscessammmneues sucuusracasunanasnnscenneensseouccacessanr) 


end, 
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{BR- 16,8 We Ned 
Program Vis_CCO; 


vVIsS_cCD 1 OF § 


uses 

ert, 

cos, 
printer, 
tpextmen, 
eedutils, 
cedrange, 
codmetrc; 


pointer; 
atringl4); 
image_neme; 
image_erray; 
integer; 
integer; 
integer; 
Integer; 
integer; 
integer; 
real; 
comment_info; 
fetop comment_info; 
expt ime reai; 
filter comment_‘nfa; 
etert hugeptr; 
entropy real; 
display video_erray; 
atd_dev real; 
ekewnens real; 
kurtosis real; 
upbound integer; 
loweound integer; 
windows window_errey; 
clutter real; 
reynoid real; 
histo histoptr; 
cont _timg video _arvey; 
mincont 
maxcont 
cont0S : Integer; 
cont 25 integer; 
cont50 : integer; 
cont?s intege:, 
ones integer; 
Junk word; 
procmonth word; 
procday word; 
prechour word; 
proc#in word; 
azimuth teal; 
elevation real; 
uv : GLVptr; 
smalligly real; 
lerge_glv : real: 
sites _fite text; 
ingtile 1 text; 
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cutfile + text; 
inge_filenama : etringt2; VI8_ccD 2 OF 5 
matr filename : string!2; 
data_loc + etringt6; 
etr2 + eteingl2); 
ent t Integer; 
tep_neme t etring(12); 
code s integer; 
view {mag 1 Integer; 
view t Integer; 
azel_tile 1 text; 

axrdeg 1 Integer; 
arin + Integer; 
arsec + Integer; 
eldeg : integer; 
elaln t Integer; 
eleec t integer; 

flag 1 byte; 

i] 1 byte; 
exp_time a tent; 


begin 
eesion(sites ¢{le, "ALL_CCO.IXT’); 
reset(sites_filed; 
readin(sitec_tile, (mge_¢ilename, etr2, metr_filensme); 
wnilec (ege_tilename <> ¢' ) do 
beatin 
ession¢ingfile, 1MGS_PILEWAME); 
reset( tmgfile): 
ceadincimefile, data_loc); 
ent 18 0; 
reading imgfile, name); 
whilet nam <> '' >) do 
begin 
ouslencouttile, METR FILENAME); 
{8I-) 
append(outfile); 
(Si) 
{f ClOreault «> 0) then rewritecoutfile); 


ent is ent * 13 

writen; 

uriteln(output, ‘Processing image #, cntr4, '  ', data_loc, name); 
GetDetec junk, procmonth, procday, junk); 

GetT imeCprochour,procmin, Junk, Junk); 


{f copy(neme,1,2) = "EG! then 
asefentazel_file, ‘€G_AZEL.1xT!) 
else 
'¢ copy(neme,1,2) ¢ ‘YG! then 
suntonCazel_ file, ‘YG_AZEL. Txt’) 
else 
If copy(name,t,2) * "AK! then 
osotgncage! file, -AH_AZEL. TXT!) 
else 
{? copy(name,1,2) © ‘AQ’ then 
eesipn(asel_file, ‘AG_AZEL. TXT") 
else 
1f ( copy(iveme,1,2) © *#0') end 
{ (copy(name, 5,2) © °26') OR (copy(neme,3,2) © '27°)) then 


 atnaEEEnnene angen 
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assigncazel_file, ‘FO_AZELI.TxXT’) 
else 
1f ¢ copy(neme,1,2) © ' FD‘) and 
( Ceopy(name,3,2) = 41") OR (copy(name, 3,2) = 12/9) then 
senigntezet_file, ‘PDAZEL2.tx1!) 
alee 
(f copy(name, 1,2) @ ‘HL! then 
essign(azel file, ‘HLAZEL.1x?’); 


VIS_CCD 3 OF 5 


TMT WE OMS Cw Linchonts ) 3 
resetiazel_fite); 
repent 
val(copy(name,7,2), viewing, code); 
readin¢azel_file, view, axdeg, azmin, azsec, eldeg, elmin, elsec) 
until ( view «© view img ); 
closecazel_ file); 
azimuth s= ardey + ((Cazsec / 60) + azmin) / 60); 
elevation 19 eldeg + (({eisec / 60) + elmind / 60); 


READ_RANCE (data loc, name); 
1¢ © (opy(neme, 1,2) © ‘YG") or (copytname,1,2) © ‘yg’? ) then 
FINO, BOUNDS_YPG(aximuth, elevation, lowbound, upbound, name) 
else 
FINO_BOUNDS( lowbound, upbound, name); 


{if ¢(lombound: upbound)<60) then 
begin 
writein¢‘bounds restriction’); 
WALT; 
end; 


WAAK(p) 3 
new (GY); 


MARKEXTMEM( Start); 

READ_IMAGE(image, name, GV, min, max, percdS, percOS, median, mode 
, @ean, room, fatop, exptine, filter, atart, entropy 
» lowbound, upbound, date loc); 


email gly re 0; 
lurge_gty re 0; 

for 1 r= 0 to 10 do 
begin 


omeli_giv c= small_gly + gv*.bins(t); 
farge_elv :© large giv + gv*.bins(16299- i) 
ond} 
emai l_olv re Comalllgiv / gv*.bime[16300)) * 100; 
largegiv c» (lergegly / gv*.bine[16300)) * 100; 


1f ¢ (Cemall_glv+large_giv) < 0.50 ) THEN 
begin 
CREATE_OISPLAY(imege, display, std dev, min, perclS, mean, skewness 
, kurtesis, lowbound, updound); 


nen(histo); 


CONTRASTA(hiato, lowbound, upbourx, image, windows, cont img 
» Maene); 
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aincont re PERCENTILE (histo, 2.0/int(histo*.bins (6192) >); 


maxcont 2 PERCENTILECHiste, 0.9995); 
Contos r= PERCENTILECHIsto,0.05); 
contzS 18 PERCHUTILECHIsto,0.25)) 
cents 3 PERCENTILECKI ato,0.50); 
Cont7S sm PEACENTILECHIsto,0.75); 
eoncOS «re PERCERTILECHiato,0,95); 


1? (trunccexptiometl0) «> 1) THEN 
begin 
enalencexp time, EXP TIME.1x1'); 
{ti-) 
eppend(exp_t (me); 
cat+) 
4f ChOresule <> 0) then rewritecenp time): 
writelnCenp time, mame, ' ', axptimerés2): 
close(exp_tima); 


CWOTE: GLY's are Normalized when the image is reed.) 


end; 
flag 0 0 


wean 38 9909p; 
ain := 9999; 
max 38 999; 
peredS = 9909; 
perc9s 15 9999; 
median 10 9999; 
@ntropy 22 9999; 
Gkewness 88 9999; 
kurtonia 1© 9999; 
mincont ¢= 9999; 
MAXCONt t= OPV; 
contdsS 12 9999; 
cont2S 18 9999; 
contsO t= 99T9; 
cont?S 12 9999; 
contes r= 9999; 
flag 19 1; 

end; 

REL EASEEXTMEM( Start); 

RELEASE(p); 


writeincoutfile, neme:0, ° ', procday:4, ' ', procmonth:3, ' 
« Procmouris, ' *, procminss, © ', azt@uth: 10:6, °° 
» whewation: 1014, ¢ ', upbound:5, ¢ ', lewboundsS, ' ¢ 


+ Wodar5); 


writelmn(outiile, means7:9, ¢ ¢, mimé, 6 ¢, pareds:7, ' 
« mediant?, 6 ', perc0S:7, ' *, max:?, | ', etd_devi9:2 


~« ' ', @ntropy:1013); 


writeln(outtile, skewwse:1013, 9 ', kurtosia:10:3, mincont:7 
oe ff, coneesr?, § *, comt2ss7, * 4, contsa:7, * & 


2 Cone7Si7, § 4, cantO5:7, ' /, mexcont:?); 
writeln(qutfile, smeli_givi2é:2, 4 ', large_glvié:2, ' 
, flagged, © 4, expt tmer6:25; 


closecouttile); 
teadiniimatile, name) 


VIS_CCD 4 OF 5 
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> 


end; 
closecimgfile); . VIS CCD 5 OF § 
readin¢aites_file, imgs_filename, str2, metr filename) 
end; 
closaceites_file) 
[NO. 
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Unit ccpkanee: 
‘ CCDRANGE 1 OF 5 


iatertace 


uses 
printer, 
codutils; 


conet 
hpin = 0,918; 
vein = 0.198; 
trarow = 38; 
mnecol # 37; 
main «© $0; 

terget_width = 6; 


type 
ronge_image © erray{!..rngrom, V..rracol) of integer; 
window.rec © record 
rarge : integer; 
hetze s integer; 
velze : integer; 
ond; 
window_erray © array(0,.tamin) of window_rec; 


ver 


ranges ; range_image; 

[ CO ee eed) 
procedure resd_range( dateloc =: stringté; 

imagename : lmage name); 


Procedure find bounds(var lowbound : integer; 
var upbound oo: integer: 
imagename : {mage_name); 


procedure find bounds _ypg( ar + reel; 
el + teal; 
war lowbound : integer; 
var upbound 1 integer; 
imagename: image_name); 


function get_range(row t Integer; 
col + Integer; 
imagenamn : Image name) : integer; 


Procedure init_windows(var windows : windew_errey); 


Procedure get_window( range integer; 


1 
windows : window_erray; 
var lowrangs : Integer; 
ver uprenge | integer; 
1 


wi rdow_ tec); 


{epi ementat ion 
(ean 


Procedure reed_range( dareloc + atring!6; 
(emgename ; iaage name); 


ver 


§, J] : integer; 


ASO 
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intitle 1 ffle ef integer; . 
junk 1 Integer; 
etring2é 2 string (26); 


begin 
{¢ vopy(imagename, 1,2) © 'EG’ then 
atring2é 1= detaloe * /EGL’ + copy(imagename,7,2) 4 ' eng! 
else 
4 copy(imagename,1,2) © ‘YG’ then 
etring26 1= dateloc + *YPG! + copy(imagoname,7,2) * ‘.rny! 
else 
{f copy( imagereme,1,2) = ‘AH’ then 
atring2é 1" dataloc + ‘APH! + copy(imagename,7,2) + '.rng! 
else 
if copy(imeagename,1,2) © ‘AG! then 
atring2é = dateloc + ‘APG! ¢ copy(Imagename,7,2) * ‘rng! 
elee 
if € copy({megename,1,2) © 'FO') and F 
€ (copy( imagename, 3,2) © '26°) OR (copy(imagename,3,2) = /27')) then 
string2s ra dataloc + ‘FID’ « copy(imegename,7,2) + 41) + ‘y rng’ 
else 
If ( copy( imagename,1,2) * ‘?0°) and 
¢ (copy  imegename,3,2) = 411") OR (copy( imagename,3,2) © 12/9) then 
String2é s© dataloc * ‘FID! + copylimagename,7,2) + /2! + ‘.eng! 
else 
if copy(imagename,1,2) " ‘HL’ then 
atring26 t= dataloc + ‘HTL’ + copy(imagename,7,2) + ‘.rng’; 


essigntinfile, strings); 
reset(intile); 
tor { 18 1 to 76 do 
for ] 18 1 to rngcol do read{infile, junk); 
for { te rngrow downto 1 do 
for } t# 1 to prgcol do read(inttle, rangeati, 3))7 
elowacintiled; 


GOES HAAN EEee OU ERAS ESE eERC ED UELA REA BNRAeeaeELeEAsEseES) 
function get_renga(cow 1 integer; 

col 1 Integer; 

imagename 1 image_mame) ; integer; 


{CONTAINS CHEAP FIM TO FILTER QUT BAD RANCE DATA) 
ver 
hewcel + Integer; 
test oo: integer; 
begin 
hewcol r= round(eol / '9.0); 
{f (rmwcol <1) then 
neweol 1 1 
else 14 (mewcol » $7) then 
mewcol 16 $7; 


wet_ronge 1° renges (row, newcol); 
test 18 ranges(rou, oeucol); 


6 ¢ ¢ (copy imeguname,1,2) © 'YG’) of (capy( imagename,1,2) = ‘yg’) ) 
end (test « 700) ) then ; 
wet_ramgere2000;  (<---FIK TO FILTER OUT SAD RANGE DATA FCR YPC) 
end; 


(Mee veensessnVeEEE TA STESUUSSSEES SNE STERROSES Cane e AM eansoneE ATO S RABE REE es) 
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. Procedure find_bounds(var lowbound 1 integer; 
var upbourd =: Integer; 
imagename 1 |wage_name); 


conat 
lower # 300; 
unper = 29000; 


var 
{,j 2 Integer; 
done : boolean; 


begin 
done 16 FALSE; 
Lie jd oie 70; 
while mot(done) do 
begin 
{f GET_RANGE({ , j, imegename? < upper then 
begin 
upbound ra {; 
done t@ true; 
ordi 
fist 4; 
end; 


Gone 12 FALSE; 
1:8 384; Joie 70; 
while net(dene) do 
begin 
if GRY_RANGEC! , J, imageneme) » lower then 
begin 
Vonbound 38 1; 
Gone := true; 
and; 
fief + 
end; 


{f upbound < 47 then upbound 1© 47; 
{¢ lowbound >» 335 then lowbound 1 335; 


{nctovus 
Procedure find_bounds_ypg( a real; 
el real; 

var lowbournd : integer; 

var Upbound in.eger; 


\magenane: imege name); 


CCDRANGE 3 OF 5 


BUeezeaanssestranunsavedtis SegueasauauE) 


(determines range limite of processed imagery to be applied 


for YUMA DEM/Val) 
const 
tower # 300; 
upper = 29000; 
type 
date pair © record 
ant integer; 
els raat; 
end; 
ver 
i,} + Integer; 
upper_done 1 buolean; 
ridge + array[1..31) of data_pair; 
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bottom 


begin 
ridgettl.ez 
ridge (2) .0z 
ridge (3) .az 
ridge (4) 02 
ridge(S}.az 
ridge (6) az 
eidget7) az 
edge (8) az 
ridge (9) a2 
ridge(10).az 
ridge 11) az 
ridge (12) az 
ridge (t3) ex 
ridge (14) .az 
ridge (15). ex 
ridgel16).az 
ridge lt7) az 
ridge ll8) ,¢% 
ridge (19) .a2 
ridge (20) a2 
ridge(21) az 
ridge (22) .az 
ridge (23). az 
ridget{24).aa 
ridgef23],a2 
ridge (26) .az 
rtege(27],02 
ridge (ze) .a2 
ridge l29) a2 
ridge (30) ax 
eidge(3{].a2 


upbound :2 1 


ce 


tied; doa 


begin 


end; 


(unvense 


var 
fs int 
begin 


t real; 


10202 
90204 
18206 
19208 
sa210 
e202 
192d 
rege 
182d 


upper done t= false; 


= 268; 


while mot (upper done) 


eger; 


ridgett).el 

ridge(2).el 

ridge(3).el 

ridgn(4).al 

ridgetS).al 

cidge(6).el 

eidgel7).el 

ridge(@).et 

ridgel?).el 

ridge (10) .el 
eidgetti).el 
cidgelt2).el 
vidge(13).e( 
ridge(14}.el 
ridue(15).¢4 
ridge({6}.e\ 
ridge(i7).e| 
eidge(16).e1 
ridge lt92.e\ 
eidge (20) .e4 
eidge(2t).et 
eidge(22).@l 
eidge (23) .el 
ridge (24).e\ 
ridge (25).el 
ridge (26) .el 
ridge(27).el 
ridge(2@).el 
ridge l29).e1 
eidge (30) ol 
eidgel3i).e4 


do 


procedure init _windowa(var windows : 


18 92.3; 
t= 92.6; 
28 92.7; 
rm 92.6; 
18 92.6; 
12 92.6; 
ra 92,5; 
rm 92,5; 
18 92,5; 
1092.4; 
2092.1; 
2091.8; 
1091.8; 
1891.7; 
1e9t.7y 
091.8; 
1091.9; 
1292.13 
8992.9; 
r091.9; 
1292.1; 
1092.0; 
1091.6; 
1091.9) 
1991.7; 
1991.8; 
1091.3; 
1091.0; 
etd; 
1097.0; 
1997.0) 


while ridgeti) axcroundiaz) do {:i41; 
bottom re cidgeti).el; 


{f GET_RANGE({, J, {magename)<upper then 
begin 


lowbound 1m round((bottom-e()#¢1.0/¢¢vpin/1000.039S7.3))) + 192; 
{f upboundes? then upoound 18 47; 
1 lowbound? 335 then lowbound: #335; 


BERENS SUSadeseAUNG eae Recor eEenUSeunennsausacseaueuas) 


window array); 
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CCDRANGE 4 OF § 


AS3 


—- +. 


unit CcOMetre; CCDMETRC 1 OF 13 


interface 


uses 

ert, 
tpext men, 
ecdutils, 
ccdrange, 
eedtgtin; 


const 
histosize = 16383; 


type 
hintoptr = *histogram; 
histogram * record 
bins : array{-8100..8192) of longint; 
end; 


t eee ere errr eet t rt irr errr r rr errr es saeeees) 


Procedure contrest(var histo s histeptr; 
lowbound : integer; 
upbdound =: integer; 
image : bmage_array; 


windows : window array; 
var cont_img : video array; 

tirat_ tpt t terget_ptr; 

imagename ; image name); 


procedure contrast2(var 


histo histoptr; 
lowbound : integer; 
upbound =: Integer; 
{nage : tmage_array; 
windows =: window srray; 
var cont_img : video array; 
\magename : image name); 


Procedure contrest3(var histo t his'optr; 
lowhound ; integer; 
upbeund oo: integer; 
image : image erray; 
windows =: window srrey; 

var cont_img : video array; 
imayename : image neme); 


Procedure contrestd(var histo histaptr; 
lowbound : integer; 
upoound =: integer: 
image : image array; 
windows window _erray; 
var cont_img video array; 


imagename : image name); 


tunctlon GTPC 


value =: real; 


histo i Bistopte) : reel; 


function PERCENTILE( war hinto : histoptr; 
percent : real): integer; 
function GLY_G1Pidark > integer; 
hist : GiWptr) : real; 
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procedure update targets(var target windows : tgt_erray; 


var ram_tots Integer; CCDMETRC 4 OF 13 
window 2 window_rec; 
Tow + integer); 


var 
i, J: integer; 
topo: integer; 
begin 
top := row ~ window. veize; 
for { :@ 1 to mum_tgts do 
begin 
1f Ctarget_windows(1},71 < top) then 
begin 
for j r= 1 to mum_tgts + 1 do 
target_windows[}) © target_windows(j + 1); 
ram_tgts 5 num_tgts - 1; 


end; 
end; 
end; 
[ Corry eevee sort ewaseeeseese Se ew ecenasasaascmanscouusseee) 
ver 
1, J, k&, | 1 {nteger; 
run_tgte + Integer; 
totlist 1 tyt_list; 
target_windows ¢ tgt_erray; 
aridetze t word; 
videotize 2 word; 
videosizek 1 word; 
accumi nner t longint; 
accumouter t longint; 
ramp! xel : integer; 
meaninner 2 real; 
mheonouter 1 
CCD Contrast 2 integer; 
offset : integer; 
range + integer; 
window 1 window rec; 
lowrange 1 Integer; 
uprange 1 Integer; 
currtgt 1 tgt_liax; 
workspace t erray(t..71) of workptr; 
procede + boulean; 
top : integer; 
bottom $ integer; 
lett 1 integer; 
right 1 Integer; 
tempi : integer; 
temptop t integer; 
tempbot tom : integer; 
wt dow | ricdex, : integer; 
current tgt i Integer; 
acur, your : byte; 


weur2, yeur2 + byte; 
wcur3s, yours: byte; 


begin 
for fie t to 71 do 
neu(workepaceti)); 
for 4 ce -B192 to BI92 do 
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histo*.binali) r= 0; 
make_tot_windows(first_tgt, target windows, cum_tets); 
sort_targete(terget_windows, num_tots); 
gridaize :* sizeof(dataspace); 
ortdsize :* trunc(gridsize / 2); 
range :* g@t_range(upbound, left_junk, imagename); 
get_window(range, windows, lowrange, uprange, window); 
writecoutput, ‘Storing Now ‘); 
xeur3 15 wherex; 
your} 1° wherey; 
top 1* upbound - window.veize > 27 
bottos 15 top ¢ 70; 
for i 1© top te bottom do 
begin 

gpotoxy(xeur3, ycur3); 

writacoutput, 1); 

moveextmen(image(i}, ptrtohuge(workspace(i - top ¢ 1]), gridsize); 
end; 
writetn(output); 
k 1) upbound; 
offset t= window.valze * 3; 
tytlist 6 nil; 
writecoutput, ‘Processing Row '); 
xour 18 wherex; 
yeur 18 wherey; 
walle (Ck <* lowbound) do 
beyin 
po.oxy(xcur, your); 
writecoutput, k, ‘ offset ‘, offser, ’ Coium '); 
meur2 = wherex; 
yeur2 t# wherey; 
check target(k, window, ietliat, rum tgts, target windows); 
range reget range k, (laft_jumk + window.heize + 1), imagename); 
(f (range < lowrange) or (range > uprange) then 

get_window(range, windows, lowrange, uprange, window); 
C18 left junk ¢ window. hetze; 
while Cl <= (mumcol - right junk - window. hsize)) do 
begin 

Procede :* true; 

curetat :* tetiist; 

left te ( + window maize; 

right te | * window.hsize: 

top te kh | witnow.vetze; 

bottom :5 k ¢ window. valre; 

while Ccurrtat*.mext « nil) do 

tegin 

current tgt i currtet* num; 
write(n(output, current tgt, ¢ /, currtgt*.num); 


4f CCtorget_windows (current _tgt).x2 > left) and (target windows[current tgt).az c= cight)) then 


procede t= falae; 


if (Ctarget windows (current _tgt).at >e left) and (target_windows [currant tgt] xi ce cight)) then 


procede := false; 


CCDMETRE 


OF 


13 


{ (Ctarget_windows [current tyt).y! >= top) and (larget_windows(current_tgt).yl «2 bottom) then 


procede := false; 


if (Ctarget windows lcurrent tgt).y2 >« top) and (larget_wiidowslcurrent tyt].y2 <= bettom)) then 


procede := false; 
currtgt te currtgt* nent 
ana; 
{f procede then 
begin 
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accumouter 1* 0; 

mumpixel :* 0; 

range :© get ranmetk, 1, imagename); 

if Grange < lowrange) or (range > uprange) then 


————__—_—_—_—_ 


CCDMETRC 6 OF 13 


fet_window(range, windows, lowrange, uprange, window): 


gotory(xcur2, youre); 
cireol; 


write(output, |, ‘ HSIZE ', window.hsize, ‘ VSIZE ", window.veize, ‘ { ‘); 


xeurd ss wherex; 
yeur3 := wherey; 
if Coffeet - window.veize < 1) then 
bettn 
tempbottom :s k - window.vsize - 1; 
temptop 18 k + window. vsize; 
for tampl :© tempoottom to temptop do 
begin 
gotoxy{xeur3, yeur3); 
windewlndex := tempi - tempbottom + 1; 
write(output, tempi, ‘ ', windowindex); 


@oveextmem |mage (tempi), ptrtohuge( workspace (windowlneex)), gridsize); 


end; 
offeet r= window.valze + 1; 
wd; 


for i se Coffeet + window.vsize) to (offset + window.vsize) do 


begin 
gotoxy(xcur3, ycur3); 
clreol; 
uriteceutput, 1); 
for j :a left to right do 
begin 


accumouter 1¢ sceumouter ¢ workspace(i)*.data(]); 


Nuapinel :s numpixel + 1; 
end; 
nd} 
eccuminner 1% 0; 


for | te Coffaet + trune(window.veize / 2)) to (offset + trunc(window.veize / 2)) do 
for j i= (1 - trume(windowsheigze / 2)) to (1 + trunc(window.heize / 2)) do 
accuminner te accuminner + workspace{()*.datelj); 


accumouter 28 accumouter = accumimmer; 


numpinel r= numpixel - (window.veize * window heise); 
maninner :© accunimmer / (window, waize * window.hsize); 


mesnouter 38 accumouter / numpixel; 
CCO_contrast :° round(meanocuter - mesninner); 


hieto*.bins(SCO_contrast) := histo*.bins (CCD contrast} + 1; 


histo*. bing (8192) := histo*. bing (6192) + 1; 


end; 

(rele? 
en; 
if Cofteet ss 36) then 
begin 


for ) ¢@ 1 to 70 do 
workspace[j) :* workspace() ¢ 1); 


@oveentmern( image(i], ptrtohuge(workapace(71)), gridstze); 


end 
alse 
offset se offset + 1; 
mie ks 2; 
updete_targets( target _winsows, mum tgts, window, k); 
end; 
ur ltelnCoutput); 
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end; 


(SeseusncecancesmaseSesenssrseenases 


Procedure contraste(var hiate 


const 
skip © 2; 


type 


ver 


lowbound 
Upbound =: 
image 
windows =: 
contling 1 
imagename : 


tat_window © record 
mit integer; 
yl: integer; 
u23 integer; 
y2.: integer; 


end; 


SRE SS ee CSREES Rees cere ERE RSEEEEEEEE Ss) 


: integer; 


integer; 


: image_arrey; 


window_array; 
video array; 
image name); 


tet_array © array(1..20) of tgt_window; 

workptr = “dataspace; 

dataspece « record 
cate : erray(1..1026) of Integer; 


ond; 


tot Uist = “tet liet_rec; 
tet_Uist_rec = record 
ram =: integer; 
next 3 tet_list; 


ond; 

ver 
4, jy kl 
fum_tyte 
oridsize 
videosize 
videosizek 
accuninner 
eccimouter 
niuarpi nel 
meaninner 
meanouter 
COO _Contraut 
offset 
range 
Wi rho 
\owrange 
Upr ange 
curetyt 
workapece 
Proc ede 
top 
bottom 
left 
right 
temp! 
teaptop 
teepbot tom 
window | ndex 
current_tgt 
aeur, yeur 


: integer; 


iaveger; 


+ word; 


word; 

word; 

longint; 
longint; 
integer; 
real; 

real; 

integer; 
Integer; 
Integer; 
window rec; 
integer; 
integer; 
tet tint 


errey[l..91) of workptr; 


integer; 
integer; 
integer; 
integer; 
integer; 
byte; 


cee ty | 


CCDMETROE 


7 OF 13 
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xeur2, yeur2 : byte; 
xeurJ, yours: hyte; CCDMETRCO 8 OF 13 


begin 
for { :2 1 to 91 do mew(workspece(i)); 
for i := -B8100 to 8192 do histo*-bins[{} := G; 
etidsize := slzeof(dataspsce); 
eridsize :* trunc(gridaize / 2); 
Tange := get_renge(updound, left_junk, imagename); 
get_window(range, windows, lowrange, uprange, window); 
xeurS := wherex; 
yeur3 := wherey; 
top 1s upbound : window.valze - 2; 
bottom :2 top ¢ W; 
for {1 1® top ta bottom do moveentmem( lmuge(!), ptrtohuge(workspaceli - top * 1}), gridsize); 
ko o:6 upbound: 
offset :© window. vefze ¢ 3; 
writeCoutput, ‘Processing Row ’); 
ucur 1© wherex; 
yeur t= wherey; 
While Ck <= lowbound) do 
begin 
@otoxy(ucur, yeur); 
write(outpuc, &); 
renge reget_rangeck, Cleft junk + window.haize + 1), (magename); 
{f Crange < Lowrange) or Cranmge > uprange) then get_window(range, windows, lowrange, uprange, window); 
(orm (ett junk + window. Apize; 
while (l= (mumcol + right junk - window.hsize)) do 
begin 
left rm + window.nelze; 
right to | + window. Wsize; 
top 12 k - wintiow.walze; 
bottom ta kh + window. vaizes 
accumouter i# 0; 
numpinet :8 0; 
range re get rangeck, |, imagensme); 
(f (range « lowrenge) or Crange > uprenge) then get_windowirange, windows, lowrange, uprange, window); 
if (offeet - window.vaize <1) then 
begin 
tempbottom ie k - window,velze - 1; 
temptop tek ¢ window. vsize; 
tae vemmi i® TeMMOTTON te ‘emtop do Moveextmam( image (temai), 
ptrtchuge(workspacetwindowi..dex)), gridsize); 
offeet :s window.vsize + 1; 
end; 
for i ra (offset - window, veize) to (offset * window.valae) uo 
for j := left to right do 
bepin 
Accvmouter :* eccumouter * workspecel!)*.deta(j]; 
Muspinel := mumpixel © 1; 
ord; 
eccuminner :© 0; 
for { 10 (offmet - trunc(window.veize / 2)) to (offeet + trunc(window.velze / 2)) do 
for } :8 Cl + trumc(window. natze / 22) to (l + trunctwindew. helse / 2)) do 
accuminner :* eccuminnmer * workapace(i}*.datal)); 
@ccumouter 15 accumouter - sccusimner; 
mampixel 1" rampixel + (window. velze * window. heize); 
@eaninner is eccuminner / (window.vsize * window. hsize); 
@eanouter eccumouter / mumpinel; 
CCO_contrast 1° roundimesnouter - meaninner); 
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{1 (CCD_contrast>8191) then Ccl_contrast:=8191; 
if (CCD_eontrast<-8100) then CCO_contrast:«-8100; 
histo*.bine(CCO_contrast) := histo*. bins (CCO contrast) + 1; 
histo*.bins{8192) := nlete*.bins(a192) + 1; 
(sa l+ skip; 

end; 

{f (offset >* 46) 
then 


CCDMETRC 


begin 
for J := 1 to 89 do workspace(j) :* workspace{j ¢ akip); 
moveextmen( image(ke44), ptrtohuge(workspace[90)), gridsize): 
moveextmem( image(ke45], ptrtohuge(workspace(91}), gridsize); 
end 
elee offaes t= offset + akip; 
k ce k + skip: 


procedure contrast3(var histo : histoptr; 
luwbound o; integer; 
upbound + integer; 
image : \mage_arrey; 


windows +: window_arrey; 
var cont_iImg : video array; 
imagename : image_name); 


conat 
skip = 2; 

ver 
i, J, 1b, Jb 4 Integer; 
rum_tyts 1 integer; 
gridusize ! word; 
videonize t word; 
videosizek word: 
accuminner + longint; 
acciumouter t tongint; 
num {net : integer; 
mean nner : real; 
meanouter : real; 
CcD_Concrast ; integer; 
range i integer; 
uindow + window rec; 
lowrange : integer; 
uprange : integer; 
top : Wnteger; 
bottom : integer; 
lett : Integer; 
right : integer; 
temp! 1 integer; 
temptop : integer; 
tempbot { om : integer; 
Window | ne x : Integer; 
xcur, yeur : byte; 
workrow : work row, 

begin 


for § re -6100 tr 8192 do higto*.bins{i} :* 0; 
Oridsise :* sizect(workrow); 

uridsiae := vrunc(gridulse / 2); 

range :* get_range(upbound, (eft junk, imagename); 
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OF 


eavasesevuaceusssean) 


13 
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get_window{range, winders, lowrenge, uprange, window); 
writeCoutput, ‘Proceasing Mow /); CCDMETRC 10 OF 13 
mcur i wherex; 
your :# Wherey; 
ire us ound; 
while (i <= lowbound) do 
begin 
gotoxy(zeur, your); 
writefoutput, 1); 
renge reget rangeci, Cleft junk © window.nsize + 1), imagename); 
If (range < lowrange) or (range > Uprenge) then get_window(range, windows, lowrenge, uprange, Window); 
J am lefty junk + window, haize: 
while (j <* (numeol - right Junk - window.hsize)) do 
begin 
range :© get_rangeti,j, imagename); 
if (range < lowrenge) or (renge > uprange) then get_window(range, windows, lowrange, uprange, window); 
left ss J + window. hsize; 
eight :" |] + window. heize; 
top is i: window. vaize; 
bottom 22 {| © window.veize; 
eccumouter re 0; 
rampixel r= 0; 
tor ib t= top to bottom do 
begin 
woveextmem( image (ib) ptr tohuge(addr(workrow)),gridsize; 
for jb :¢ left to right de 
begin 
eccumouter :* accumouter + workrowljb); 
rampixel 25s mampixel + 1; 
end; 


and; 
eccuminner := 0; 
for ib :® (i + trunc(window.vsire / 2)) to (ie trunc(window.velze / 2)) do 

begin 

moveextmem( image(ib) ptrtohuge(eddr(workrow)),gridsize); 
tor jo se () + trunc(window.nsize / 2)) to (J + trunc(window.heize / 2)) do 
eccumiimer :* accuminner + workrouljb); 

end; 
accumouter :* accumouter + accumimer; 
mempixel := mampixel + (window.vilze * window, neize); 
meaninner := accuminner / (winden.veize * window.hsize); 
@eanouter := arcumouter / rumpizel; 
CCO contrast :* round(mennouter + meeninner); 
if (SCD_contrast>8191) then CCO_contrest:28191; 
if (CCD_contrase<+8100) then CCO_contrast:=-8100; 
histo*.binsiCCO_contrest) 1= histo”.bins{CCD contrast) + 1; 
hiato*.bina (8192) += histo*.bing (8192) + 1; 
Jo:e j ¢ akip;  (colum skip) 

end; 
fore i + ekip; (row akip) 
ard; 
end; 

CUCU REC ces TENET REE SRE CR ERSTE RE AAU RNGRATET OSES LSRNS SRRSCEE REESE TSETETER ERE RAE EEE AReAe Fee Has ER INTAKE EERIE ATA) 
procedure comtrasta(ver histo : histoptr; 
lowbound oo: integer; 
upbound : integer; 
image : lange array; 
windows : window array; 

ver cont_imp : video_array; 

imagename : image name); 
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const 
skip 3; CCDMETRC 11 OF 


ver 

i, j, 1b, Jb: integer; 
Pum_tgts : integer; 
aridsize word; 
videosize : word; 
videos lzek : word; 
accunimmer + longint; 
accumuter : longint;: 
rumpinel : integer; 
inmerpixe(s : integer; 
outerpixets iateger; 
mroninner teal; 
meanouter : real; 
CCD_Contrast integer; 
range : integer; 
window : window cec; 
lowrange integer; 
uprenge (Integer; 
top integer; 
bot tos integer; 
left 2 integer; 
right integer; 
tempi : integer; 
temptop : integer; 
tempbot tom : Integer; 
tebindow inde x : integer; 
acur, yeur byte; 
workrow work row; 


begin 
for | :@ +8100 to 8192 do histo*.bins[i] :* 0; 
gridsize © aizeat(workrow); 
gridsize :® trunc(gridvize / 2); 
fenge :* get_range(upbound, left Junk, Imagename); 
get_window(range, windows, lowrange, uprange, window); 
write(output, ‘Processing Row ‘); 
xcur t= wherex; 
ycur := wherey; 
1 :@ upbeund; 
while (i <* lowbound) do 
begin 
gotony(xcur, ycur); 
write(output, 17; 
renge :sget rangeti, (left_junk ¢ window.haize ¢ 1), imagensme); 
{¢ (range < (owrange) or (range > uprange) then get_window(range, windows, lowrenge, uprange, window); 
i re left junk ¢ window. Weize; 
while (j <= (mumcol - right_jumk - window.hsize)) do 
begin 
range t© get_range(i,), imagename): 
if (range < lowrenge) or (range > uprange) then get_window (range, windows, lowrange, uprange, window); 
feft t# }| - window. hsize; 
right :© ) * window.haize; 
top := 1 - window. valize; 
borrom :8 4 ¢ window.vsize; 
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ee ee do ; CCMMETRC 12 OF 13 


moveextmam( (mage [ib] ,ptrtohuge(addr(workrow)),gridsize); 
jo := left; 
while (/b<eright) do 
begin 
accumouter t= accumouter + workrow{Jb); 
tumpixel :* tumpixel ¢ 1; 
jbts jb + skip; 
end; 
(ore ib + skip; 
end; 
eccuminner :* 0; 
innerplxeies= 0; 
ib ra | > trunc(window.vaize/2); 
while Cib<e(i+ trunc(window.vsize / 2))) do 
begin 
moveextmem( image (ib) ,ptrtohuge(addriworkrow)), gridsize); 
Jo 18 j = frunc(window.heize/2}: 
while (Jb«m(] + trunctwindow.helze / 2))) do 
begin 
accusinner re sccuminner ¢ workrow(jb); 
innetpixelez= imerpixeis + 1; 
Jb := Jb + skip; 
end; 
ib :2 1b + akip: 
ond; 
eccumouter :» sccumouter * accuninner; 
outerpixels:= numpixel - innerpinels; 
meeninner s= accuminner / Inmerpineis; 
@tenouter := accumauter / outerpixe(s; 
CCO_contrast :§ round(meanouter > meaninner); 
if (CCO_contrast>8191) then CCO_contrast:=8191; 
{f (CCO_contrast<-8100) then CCD_contrast:=+8100; 
nisto*.bina(CCO_contrast) :* Wisto*,HinsiCCO_contrass) + 1; 
histo” .bins{8ive] re histo*.bins{6192) + 1; 
J 2 + akip; (column akin) 
end; 
fre f+ ekip; Crow skip) 
end; 
end; 
(ARETE T SOE O SRN TO SOON NEs Ln NESE ANSS SETHE NESE NES SO DNs a ONSHORE SCENE SCECELPSSR SORES A USUN UE PAU SEE EERERE ESSE CENA UREtaReesanes) 
function GIP¢ value : real; 
histo t hietoeptr) : real; 
(for contrast histogram only) 
ver 
eceum =: longint; 
temp: real; 
i 3 integer; 


begin 
eccum :* 0; 
for i :@ 8100 to round(vaiue) do accum rs accum ¢ histo*.bins({ll; 
GTP :# accum / histo*.bins (8192); 


Sversueuanaganvesense) 
fumnetion PERCENTILE<var histo : hiatoptr; 
percent : real): integer; 
var 


Long (rit; 


stop_et: 
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i 2 integer; 
accum oo: to agint: 
begin 
stop_at t= round(percent*int(histo*.bins (8192) }); 
eccun ;8 0: 
f s °8100; 
while Caccum<stop_at) do 
begin 
eccum 1=@ accum ¢ histo*.bing(t); 
friet: 
end; 
PERCENTILE :@ i; 
end; (percentile) 


(CoanunesspemenenseusssetonensconsasnsrusecenespseeuanseEssnstvaceranaseusncannecuness) 
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tunetion GLY_GTP( dark : integer; 
hist : GLVptr) ¢ real; 
(computes GTP of darkest pixel on target) 
ver 


accu =: longint; 
tap : reel; 


i 3 integer; 
begin 
eccum 18 0; 


for 1 3 0 to dark do accum tw accun + hiet*.bins{i]; 
temp r= accum / hist. bins (16300); 
GLV_GYP := 1.0 + temp; 


ord 
COUEAECGhased nao ceneunneseuomousEnssOEMEeeUeueoent 10 seve aE nay sone kn ous PERSE ROME ER eR eee teNeereassentsnankennseyers wunmeer) 


end. 
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unit CCOCOMPT; {not working, clutter valuesratdey mhicn just cannot be) CCDCOMPT 1 OF 4 


interface 


uses 
ws, 
ert, 
eedutils, 
eedrange, 
tpextnem; 


= 16.0; 


e = 0.000118; 


procedure compute_reynold( image 


ver clutter 
var reynold 
etd_dev 


Procedure campute_reynoid2( image 


(eva m neem en ae mens wane mese ened names sams eee nemssnaasasassnass) 


imptementation 


Vosibound 
var clutter 
var reynold 

atd_dev 


t 


fmage_array; 
real; 

real: 

real); 


s imagelerray; 


integer; 
integer; 
real; 
real; 
real); 


PPP erer rrr reer etre rr treet etree reer ree rere es } 


(eexeacvanecessac cern eee Ver SRUAS SCORE SCRSOREECODESR SA TEOR TEENA SseEUEpAseRneUsAUEDORSSAO) 


Procedure compute _reynold( image 


type sareeal_erray 
type sqrint_erray 


var size pix 
width 
maxrow 
mancol 
residual 
rowekip 
col_ekip 
fow_sqr 
¢ol_aqr 
fiest_row 
lest row 
tiret_eol 
Vast_col 
‘J 
eriduize 
sar_eum 


sqr clutter 


work 
a,b,c 
wmeoiver 
acum 


var clutter : 


ver reynold 
atd_dew 


= arrey (1.,50, 1..50) of real; 
= array (1..50, 1..50) of integer; 


real; 

integer; 
integer; 
integer; 
integer; 
integer; 
integer; 
inteper; 
integer; 
inteyer; 


: integer; 


integer; 
integer; 


+ integer; 


word; 
agrint_erray; 
sarres!_errey; 
work row; 
real; 

real; 

real; 


: integer; 


image_array; 
real; 

veal; 

teal); 
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range := get_range( 192,284); 
size_pix r= range * angle; 
width s# trunc( dim / size_pix); 
maxrow :® rumrow div width; 
residual :* numrow mod width: 
tow_skip :* residual div 2; 
@axco| © mumcol div width; 
residual :2 mumcol mod width; 
col_skip t= residual div 2; 
aridsize 1 siaeof(work); 
eridsize s= trunc(gridsize/2); 


for 1 :# 4 te 50 do 
for } = 1 to 50 do 
begin 
eqr_eum({,]) :© 0; 
ear_clutter(i,j} :2 0 
end; 


row_eqr re 4; 
eol_eqr t= 1; 
first_row t© 1 ¢ row_skip; 
Vant_row :2 width + row skip; 
repeat 
for 4 1* first_row to ast row do 
begin 
moveentnem( image(s) ptrtohuges acdr(work)),gridsize,; 
first_eol :2 1+ col_skip; 
\ast_col re width + col_skip; 
tepeat 
for j i= firat_col to lest_co. do 
Sqr_eumlrow_sqr,col_sqr) re sqr_sum{row_sqr,col_sqr] + work(J); 
firat col r= firsteot + widt?; 
(aat_col re last_col ¢ width; 
sol sar se col_sqr + 1; 
until¢ col_aqr > mancol ); 
col_sgqr se |; 
end; 
flrat row re first_row * with; 
Vast_row se leet raw + width; 
fow_eqr t© row_sqr ¢ 1; 
until rowagr » mexrow ); 


rowleqr te 4; 
col_eqr 18 1; 
first row is 1 ¢ row skip; 
(est_row :e width ¢ row skip; 
repeat 
for i © first_row to (ast_row do 
begin 
@ovaextmem( image(i) ptrtohuge(@oie(work)), grigsize); 
tiretcol se 1+ col_skip; 
lest_col te width + col_skip; 
repest 
for j te firat_col to lest_ce. do 
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byte_dat byte_carrier; 
real_dat real_carrier; 
ofre_pix : real; 

width integer; 
vert boxes integer; 
horiz boxes integer; 
residual integer; 
row_stert integer; 

col start integer; 

I,J, 1b, jb, k, ko: integer; 
eccun 2 reals 

range 1 integer; 
work + work tow; 
gridsize : word; 


range 15 get_range (192,288); (based on center image range) 

size pin t= range * angle; (nominal pixel size in meters) 

width sm» trunc¢ dim / alze_pixd; (pixel dimension of unit box) 

vert boxes rm (( owbound:upbound) div width; <mumber of boxes in vertical direction) 
fow_atart ¢@ upbound; 

horiz_bones = (mumcol-right_junk) div width; (number of boxes in horixuntal direct ton) 
col_stert :# 41; 

erideize 17 sizeotcwork); 

otldsize re trunctaridslze/2); 


(compute and store variance af each box) 
kimQ; {bun counter) 
for it=l to vert_boxes do (boxes in vertical) 


for ]:e1 to horiz_boxes do {baxes in horizontal) 
begin 
kb:80; (plael counter in box) 
(load carrier stray end compute box variance) 
for threcdi-V)*width+row start) to (i*widtherow start) do 
begin 
moveextmem( image (ib) ,ptrtonuge(eddr(work)),gridsize); 
for JbreC¢}-1)*width+col_etart) to (j*width) do 
begin 
kbiekbe1; 
byte_dat(kb) := work{}b); 
end; 
ered} 
brake; 
byte det (0) :ekb; Cpiaely in box) 
real_dat(k) t= VARIANCE(byte_dat); 
end; 
accum r= 0.0; 
real dat (0) rwineck); 
for kisi to riwnd(real_dat(0l) do accum:saccumreal dat (kh); 
clutter :# aqrtteccum/real dat i0}); 
reynold : (etd dev - clutter)/std_dev; 
end; 
( Seewen seen eeE Os CUNT OSAREADE SONS TOSURICCEEVUSER ESSER EEE RET ESTO RSCTA RRA KER) 


end. 
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VISITED AZIMUTH 


APG 37,00 
APG 39,50 
APG 42.00 
APG 44,50 
Apa 47.00 
APG 49,50 
APG 52.00 
APG 54,80 
APG 57.00 
APG $9.50 
APG 62.00 
Apa 64,50 
APG 67,00 
APG 69.50 
APG: 72.00 
APG 74.90 
APG 77.00 
APG 79,50 
APG 82,00 
APG 84,50 
APG 87.00 
apa 69.50 
APH 297.83 
APH 300,33 
APH 302,83 
APH 305.33 
APH 307.83 
APH 310,33 
APH 312.83 
APH 315.33 
APH 317.83 
APH 320,33 
APH 322.83 
Eau 95.00 
EGL 97.50 
€at 100.00 
EGL 102.50 
€Gh 105.00 
EGL 107.50 
cau 110.00 
EGL 112.50 
EGL 113.00 
EOL 117.50 
EGL 120,00 
EQ. 122.50 
€0L 125,00 


127,50 


Site (OEGREES) 


SCENARIO 


ELEVATION 
(DEGREES) 


91.17 
91.17 
91.17 
91.17 
O17 
WNA7 
91,17 
91.17 
91.17 
91.17 
917 
91697 
WN.17 
91.17 
91.17 
91.17 
VI1.I7 
WAZ 
WA? 
WA7 
W17 
W117 
90,75 
90.75 
90,75 
90.75 
90.75 
90.75 
90.75 
90,75 
90.75 
90.75 
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E2 


Pinclude <stdio.h> 
#irclude <fentl.h> 
include <etring, b> 
Winelude ®readi ine.c* 


aainarge, argv) 


int argc; 
char *arev(); 


¢ 
int ¢7 
cher varname(20}, (ine (100) ; 
floet soler ; 
float airteap : 
float reihum H 
float veget ; 
print¢(Ménter velue for solar An"); 
scant ("xt" ,deotar HH 
prinet(Minter velue for eirtemp AW: 
ecant ("X¢" ,belrtemp 
printt<#enter veiue for rethum \n"); 
scent (¥X%¢* Srethum % 
print? ("enter velue for veget An"); 
ecant(*X¢", dveget % 
TREE! 


print? C#\A\n\n\ninyayniatatninialninininintain); 


WODE 1! 

print t(MWODE @1\n"); 

1(eoier <x 1,4100000000E+02) goto HODE2; 
@lae goto NODES; 


WODE2: 

PRIAtHC*MODE #2\n"); 

if(veget «<#  8,5300000000€ +01) goto WOES; 
else goto TNODES; 


WODE 3: 

prinet(MMQde #3\n"); 

if (aolar <# 1,080U000000E +02) goto TWODE 1; 
else goto TNODE2; 


WODES: 

PrInttCMQDE B4\n"); 

(ftalrtemp «cs 2, 7200000000+01) goto MODES; 
elae goto WCOE21; 


WODES : 

printhCMNODE #\n"); 

{f(alrtemp <1 9.0000000000E +00) goto MODES; 
elae goto NODEI0; 


WODES: 

printfCMMODE 86\N"); 

if¢relhum <» 6. 3S00000900E+01) goto TMODES; 
else goto NODE7; 


CCVALLOA 


cev.c 1 OF 5 


re 
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WODE?: 

print¢c@NOoE a#7\n"); 

if(veget <= 9.6500000000E+01) gato NODES; 
else goto THOOES; 


MODES: 

print? (MnNaoEe #O\n"); 

{#(molar <= &.46000000008+02) goto TNOODES; 
else goto MODES; 


WOODED: 

print (Made #9\N*); 

{f¢rethum «a 6,65000000008+01) goto TNODES; 
elae goto THODE7; 


WOOE 10: 

PrINCICMNODE #10\N"); 

If¢relhum <s &,050000000NE+01) gotc NOOE11; 
else goto MODEI2; 


WODE 44: 

PrinttCMMODE #11\n"); 

(¢¢relhum <= S.SSOOQQ0000E+D1) goto THODED; 
alse goto THODEIO; 


NODE 12: 

PHINCHCMMODE @IZ\N"); 

{f(veyet <= 7.96000000006 +01) goto NODES; 
@lee@ goto NNDE1S; 


WODE TS: 

PRIMCICHNODE BT3\n#); 

(f¢veget «# 7, 7600000000€+01) goto NODEI4; 
else goto MODEIS; 


WODE TA: 

printecMnde #lain"); 

{f{alrtemp «= 2,6800000000E+01) goto THODENT; 
else goto THODEI2; 


WODE 15: 

Primt#C@MQDE MIS\n"); 

itCairtemp <e 2,$200000000€+01) goto TNODENS; 
@lae goto THODEI4; 


WODE 16: 

PRIAtPCPMODE #16\n"); 

Hecveget <e 9, 7O00000000E +01) goto WODEI7; 
else goto TmobEeZd; 


NODE ?: 

printeC*MODE #17\n"); 

(f(rethum <« &.S500000000E*01) gote TwODE15; 
else goto NODE1B; 


WOOE 18: 
Print¢(*MIDE BIB\N"); 

if(rethum <= §, TSOQ000000E¢01) goto THODET6; 
alae goto NODEND; 
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MODE 19: 


Prince ("NODE MIP\N"); 
(f(weget <# 8.3$00000000€+01) goto NADE20; 
eles goto THUDEI9; 


HODEZO: 

Prine#(*WODE #20\n"); 

if(veget <" 8. 1000000000801) goto TNODE17; 
elae goto THODETA; 


HOOEZ 1; 

PRIMEFCPNODE #21\n"); 

If (veget <= 8,0200000000E+01) goto NODE22; 
@ise goto NODE27; 


MODE22: 

PrinceCMMODE #22\r")> 

(f(eolar «= 8, 7100000000€+02) goto wane23; 
elae goto MOOE26; 


MOET: 

PrintfC*NODE €23\n"): 

Vé(relhum <=  1,7300000000+01) goto TNCOE21; 
@lee goto WOOE24; 


WODE24: 

prints ("MODE #24\n%); 

{f(veget <= 4,05000000008901) goto THOUEZ2; 
else goto WODEZS; 


HODES: 

PRIACECMNODE W25\r*); 

(f(veget <= 7.950000000!' +01) goto TMIDEZ3; 
else goto THODEZ4; 


MODERS: 

Print#CMMODE #26\n"); 

(f(veget <= 7,6300000000€+01) goto TNODE2S; 
alee goto TMODEZ6; 


WODE27: 

Printh( "MODE #27\n"); 

if(airtemp <e 3,4500000000E+01) goto NODEZS; 
else goto IMCDES2; 


WODE28: 

Prine#C™MODE #28\rr™); 

If (veget <m 8, 7800000000801) goto NODE29; 
alae goto NODES; 


WODE291 

PCINtt ("MODE #29\n"); 

(ffrelhum te 4, 2500000000E+01) guto wODE30; 
elae goto THODE29; 


WODESO: 

PrIMtt(MWODE M30\n*); 

if(veget «<= 8,55000000008~01) gota TmIDE27; 
alse goto THODEZ8; 
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COVALLOA 
CCvV.C 4 OF 5 


WODES1: 

PrlAttCMNODE #31\n"); 

{flsoler <= 1,8800000000F+02) gate TNQDESO; 
else goto TNOOES!; 


TWwODE 1: 

THODES : 

THODES 

TNODES: 

TWODES: 

TWODED + 
THODE28: 
THODE29: 
THQDE3O: 
TWODES2: 

prince c@\nin®); 
Princt(™Clees Min"); 
goto ENO; 


TWOOE2: 

TWODES: 

TWODE 6; 
TWODE20: 
TWODER' ¢ 
TWODE23 : 
TNODEZG: 
THODER?: 
TWODES 1 : 
PRIneeCM\M\AH); 
peintt("Ciass #2in"); 
goto END; 


TWOOE71 

THODE TI: 

THODE 14: 
TWODE'S | 
TWODEI7: 

THODE TO: 

TNODE 22: 

TWODE 26: 

TWODE 25; 
printi(™\n\n*); 
printh(*Class s\n"); 
goto END; 


TWODE 10: 

TWODE12: 

THODE 13: 

TNODE 18° 
Printf(Aynin"); 
Prinet(*Class Min"); 
goto ENO; 


EW: 
return); 


int (pin¢rum, (ist 
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E6 


int ram 
cher Wiett); 


¢ 
char nl20); 
char 0(20); 


stoa(rum,n, 10); 

strepy(s,*."); 

streat(s,n); 

Gtreat(e,“."); 

(f (atrlentatratr¢iiat,e)) » 0) 
return(1); 

else 
return(0); 


CCVALLOA 
cev.c § OF 5 
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finctude <stdio.h> 


Minclude <fentl.h> 
Cinclude cateing.h> 


Cinclude “readl ine.ct 


cher *argv(); 


tnt ¢; 


char varname (20), |ime(100); 


float reine 
float airtemp 
float veget 
float winddir 
float erdd 
float er30 
float ari20 
float solar 
float range 
float windepd 
float green 


Print(cSenter value 
tcant(*Sf", 8relhum 
printh(“Enter value 
Beant (*Rt", bei rtemp 
printt(“Enter value 
scant ("%4" &veget 
PrincecmEncer vatue 
scant (*R", deinddir 
printf “Enter vaiue 
scanti "Xf", £er60 
trintt(“Enter value 
acent(*St" &ar30 
BrinteCMEnter value 
econtc™ AF ber 120 
printt(“Enter value 
seanfC"Zt™ Eeolar 
printt("Enter velue 
ucant(*2/" grange 
PrinteC*Enter value 
ecent("S4", dwt dap 
Printf(“Enter value 
acent(*Nt", green 
TARE: 


Print t(AVAVANARIA\ ARRAN \n\nin\ninin\n\ntnin®); 


WORT: 
pPrintt(MNODE WI\n"); 


tor 


for 


for 


for 


for 


for 


for 


tor 


for 


for 


for 


\n"); 


An); 


An"); 


An*); 


An"); 


\n"); 


An"); 


\n"); 


\n*); 


\O"); 


\n"); 


if{veget <9 9. 3100000000E +01) goto H@EZ; 


wlee goto MODES?; 


WOOE?: 
Print¢(*WODE M2\n%); 


4#(6r0 <= 5.4 700000000E +02) goto NODES; 
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E7 


E8 


PS eg 


elne goto HODESE: 


WOOET: 

print t(PNODE @#3\n"); 

(ftqwinddic <e 2.610000G000E +02) guto WODES; 
alse goto MODE29; 


ODES: 

Printé(“MODE Bh\n"); 

{¢(rethum <= 9. 3500000000E +01) goto WODES; 
alee goto MODEZS; 


NODES: 

printtcPHoDe #510"); 

{Cwindepd <e 4. PSCQDOQOOUDE+00) goto HODES; 
elne goto WODENS; 


WODts: 

PRIMthCMMODE @6\n"); 

(f(winddlr <= 1,2500000000E+02) goto WOOE?; 
else goto MODES; 


WODE?: 

PRINCECOMODE M7 \rr*); 

{fceolar «= 5.6800000000€+02) goto MODES; 
elae goto THODE7; 


WODER: 

printt(*WOOE M&\n"); 

{f(winddie «= 4, 4500000000E+01) goto MODES; 
alee goto NODETO; 


WODED: 

printt(#NOOE @9\n), 

Ifceolar <e 1, 3100KO0000E401) goto TNODE!; 
elee goto TNODE2; 


WODE TO: 

PRIAFCMMODE M1O\N"?; 

{f(rethum <» 6,150)0000000E 001) goto WODE 11; 
@lse goto THODES; 


WOODEN: 

print#(@Mooe @11\n"); 

{f(relham <e 3. F500000000b401) gore Yat; 
ene goto NODE12; 


WODET2: 

printt(eMade @i2\n"); 

{#(veget <2 B.000Q00000E+01) goto THODES; 
@lae goto TNOOES; 


WODETS: 

printfCHMOOE @#t3\n"); 

{¢(soler «2 3. 7600000000E +01) goto WODE 14; 
else goto TNODE10; 


WODE 14; 
printt(MWOOE #14\n"); 
i¢(airtemp <e 5, OQOQ0QU0E+00) goto TNODES; 


co TALLOG 
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else goto TNHODED; 


MODE15: 

printf(*"NOOE #15\n"): 

4f(art20 <= B,GOQOOD0000E+00) gato MODE 14; 
alse goto MODES; 


MODE 16; 

printhCMWQDE W#t6\n*); 

i¢(renge <= 5. 5500000000E+02) goto MODE17; 
else goto WODE18; 


WODE17: 

pringtCMMQDE MI7\N"); 

If¢airtams <m 6. 9500000000E+00) goto THODE!!; 
e@lae goto THODET2; 


WOE 16: 

prine’(*WOOE #18\n"); 

{fCrethum <# 5.5000000000E+01) goto rMoDe13; 
e@lee goto TNE; 


WODEI9: 

Prined(MMOOE #19\n"); 

{ftuinddlr <a 2, 1300000000€ 402) goto MODE20; 
else goto MODE23; 


WODE20; 

PrinttC*HQOE #20\n"); 

1¢ (ere <e  4.1100000090E*02) goto TNODE15; 
lee goto MODE21; 


WOOE21: 

Print#C“MODE #21\n"): 

it (windepd <9 $,1000000000E+00) gote NODE22; 
else goto TNODE1B; 


WODE 22: 

printtC@MODE #22\n"); 

{t(rethum <s 4. 1500000000E+01) goto THODE1S; 
elae goto TWODEI7; 


WODE23: 

PrINtHC*NODE #23\n"); 

{fCrethum te 5, 15000000008 +01) goto THODE19; 
elme goto THOOE20; 


MODES; 

Printt(*WODE #24\n"); 

if(windapd <2 2.250000U000€+00) goto NODEZS; 
else goto TNODE2S: 


WODE2S : 

printt("WQ0E #25\n"); 

it¢walar <= S,.SS00000000E+00) goto MODEZ6; 
elue guto MODE2E; 


MODE 26: 
PriNtCMWODE #26\n"); 
ff(relhum sa 9.4$00000000€+01) goto THODE21; 
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else goto MODEZ?; 


WOCE27: 

PrINtH(SMODE #27\n"); 

{fCrelhum <* 9, 7500000000E +01) goto THODE22; 
elee goto THOOE23; 


WODE24: 

Printt(MOOE M28\n"); 

(fC winddie <@ 6, 3000000000E+01) goto TMODE24; 
elae goto THOOEZS; 


WODE29; 

PrINttC MODE M2P\N"); 

ifCulnddir <e 3, 200Q000000€+02) goto MODESO; 
else goto NODE41; 


WODESO: 

PHIATICONCOE W3O\N"); 

(fCerd0 «© 1.00000000008+00) gote woDES!; 
eine goto MODE32; 


MODES; 

print¢(MODe MIAN"); 

{t(renge <= 3, 0500000000€ 02) goto ThODed7; 
else goto THODE2B; 


WODF32: 

PrIMttCMMODE M32\n"); 

ifdranpe <e 2, 9500000000E+02) goto TNODEZ?; 
else goto NOOES3; 


WODE33: 

PRINCICHNCDE #33\n"); 

V#Cairtemp <=  1,4400000000E+01) gato WODE34; 
elee goto NODES: 


Hots: 

PrintéCMWOOE #34\n"); 

4#(ar120 <e  1,9700000000E+02) goto NQDES5; 
alae goto THODES; 


WODESS: 
PrIAth(MQDE #¥8in"); 


ifGulnddie <= 2.7100000000E +02) goto TMODESO; 


else gotu MODES; 


WODEM: 

prinktt(*CDE T34\n"); 

{f(relhums xe 9. 75000000008+01) goto NOOES?; 
eee goto TNODENS; 


WODES?: 

PrinteCeMODE H3T\N"); 

{f(ersd <= 7, 0000000000600) goto THODES!; 
lee goto MODESE: 


woDe 38: 
PrIATECPNODE #34\n"); 
{f(aedd <2 1, 0500000000E+00) gots THODES2; 
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elae goto NOOES9; 


WODES9; 

printt(*#ODE #39\n"); 

if(relhuw <n &,3500000000E+01) goto 1NCOESS; 
else gota TNODES4; 


WODE4O: 

PriMtt{#WODE #40\n"); 

({(veget <= 8,9700000000E401) goto THODES7; 
elee goto THODESS; 


WODES1; 

PrIACH(*MODE #41\n"); 

1t(windaepd <= 2, 1000000000€ #00) goto WODE42; 
elae goto MODES; 


WODES2: 

printé(*WOOE #42\n"); 

if¢(eirtemp <# 3.0000000000E+00) goto TNODES9; 
else goto THOOESD; 


WODES$: 

PrIMtfC@NODE MI3\n"); 

1{#Car120 <= &,QOQQ0O0000E+00) goto TNODES!, 
else goto MODES; 


NODES: 

Printf(MNQDE @44\n"); 

if(rethum <2 §,8S00000000E+01) gotc TRODEG2; 
alee goto HODESS; 


WODESS: 

printt(MMODE #45\n"); 

ifCrelhum <= G6. 3O0DQ00000E +0!) goto THODESS; 
elae goto NOOESS; 


WODE 44: 

PrINtf("WODE @4O\n"); 

if(windepd te 6, 3OQ0UOO000E +00) goto TCE G4; 
else goto NODES7; 


WODES?: 

printt(“WOOE @47\n"); 

if Cveget <2 &.9600000000E*01) juto TNODE4S; 
else goto TNODESS; 


WODE 68: 

printtc*MODE #46\n"); 

f(wireepd <8 2. 9000000000E +00) goto NODELD; 
else goto NODES3; 


WODE 49; 

PrINt((PMODE #49\N"); 

if(veget <= 6. 2700000000F+01) goto NODES; 
else goto THODES!; 


WODESO: 
printeCHWOOE #50\n")>; 
if(eirtemp <2 2,8700000000E +01) goto WOUEST; 
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El2 


else goto MODES2; 


WODES1: 

Print¢C*MODE #5 1\n*); 

{€(er60 t= 7, 6500000000E*02) gata THODES7; 
alse goto TRODE4S; 


WODES2: 

printé(*WODE #52\n"); 

if(airtemp <a 3.QP0Q0Q0000E+01) goto TNODESS; 
e@lee goto THODESO; 


WODES3: 

BriAtt(*WOOE #53\n"); 

if (aré0 <= 7,4300000000E+02) goto NODESS; 
lee goto TNODESS: 


WODES<: 

printeC"HODE #54\n"); 

1(er60 <m §,5200000000E +02) goto THODES2; 
alee goto ODESS; 


wopess: 

br intt("NQOE 855\n"); 

{fiwinddir «= 2, 0000000000€+02) goto NODES6; 
@lae goto THODESS; 


WODESS: 

POIMeeCMMODE M5A\N"); 

(f(range <= 2, 6500000000E+03) goto TNODES3; 
else goto THODESS: 


WODES7: 

prinet(MMOOE #S7\n"); 

{F(range «© 3, 200000000NE+02) goto NODESE; 
else goto THOOESS; 


WOOESS; 

PrInt{CMWQDE MBB\n"); 

{¢(winddir <a 4,31000000008+02) goto NODES9; 
@lee goto NODEGO; 


WOOES9: 

PrINttC™NODE M59\N"); 

(fdwinddir <a &.2500000000E+01) goto TMODES7; 
e@lee goto TNODESS; 


NODE 60: 

Prinrt(*WODE @60\n"); 

{#(er30 <= 4, 6OQ0000000E+00) gotn 'WODES9; 
wlae goto MODEG!; 


WODES1; 

princt(*WQDE @61\n")>; 

{@Cae3O <= 9, PIOCONOODDE +02) goto NODTA2; 
lee goto TNOOES3; 


WODE62: 
print((*WODE M62\n"); 
{f(aolar «#8 3, MOOD0Q000E +02) goto MODESS; 


> 
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e@lae goto TNOOES2; 
CST.C 7 OF g 
WODEA3: 
printf ("WODE W63\n"); 
if(sotar <= 1. 770000000GE+02) goto TWODE6D; 
elee goto TNODES1; 


TWODE12; 
THODE 18; 
THODE26: t 
TWODE2 1° 

TWODE23; 
THODE2S : | 
TNOOE 28: 

THODE 34: 

TWODES1: 

TWODESE: 

TWODES1: 

THODE 64: 
Priatec*\ninm); 
Printéc*Class #1\n"); 
goto ENO; 


TWODE2 : 
TWOOES: 
TWODES: 
TWODE 10: 
TWODE 14: 
THODE 14; 
TMODE 1S: 
TWODE 17: 
THODE 19; 
THODE24: 
THODE26: 
TWODE29; 
TWODESS 
THODE42: 
TMODESS: 
THODE SS; 
THODEGS : 
TWODES7: 
TWQOE GO: 
TWODE 62: 
PRINT ECR LN); 
printh(*Cl eee #2\n*); 
goto eNU; 


THODES: 
TNODE 4 : 
TWODE?: 
THODED: 
TWODE 13: 
THODE 22; 
TWODE2?: 
TWODE31; 
TWOOE 34; 
THIDES?: 
THODESO; 
THODES3: 
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THODESS : 


THODES?: 
TWODESS: 
THODES 1: 


TNODE SS: F i 
TWODESS : 
TRODESS : \ 


TWOOESS: ' 
prineéc™vn\n"); 

printf("Class @3\n"); 

goto END; 


TWOOES : 

TWODEB: 

TROOE 16: 
TwQOE30: 
THODES2: 
TWOOES3 : 
TMODE34 : 
TRODES9: 
“TWODESO: 
TWODES2: 
TMODES3: 
TWODESS: 
prinefceynine); i 
printfc"class #4\n"); | 
goto END; 


Ewo: 
return(0); 


int tsin(num, List) 
int mum; 
char (ist (); 


char n(20); 
char s{20}; 


jros(nun,n, 10); 

stecpy(s,"."): 

etecat(s,n); 

streat(s,"."); 

if (strlen(stretr(list,#)) > 0) 
return(1); 

else 
retucn(0); 


~ 
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CCTALODA 
1 OF 7 


cetdio.h 
finclude <fentl.h> 
Circlude <eteing.h> 
include "readline.c” 


ccrT.c 


main(ergc, argv) 
int argc; 
char *arevl); 


C 
int ¢; 
cher Vvarname (20), (ine(100); 
float relhum H 
float veget ; 
float sré0 : 
float solar y 
float renge i 
float #320 ; 
float winddir : 
float 630 3 
fiost sittesp 1 
float windspd i 
float green : 
printeCMEnter value for rethun Ww"): 
tcanf "Xt", 2relhum % 
printec"Enter value for veget \n*); 
scant (#Zt", dveget oH 
print¢("Enter value for eréd An"); 
scant (Xt, bsré0 Hh 
print? ("Enter velue for solar \W); 
econf. "2", solar y 
Print¢("Emter value for range An"); 
scant (*X" , Srange v 
print?(*tnter value for ar120 An); 
scant (WX, ber120 oH 
OriAte(MEnter velue for winddir \n"); 
scant (*Rt", dwinddir % 
printéc*Enter value for sr30 \n"); 
ecand( RE" bers % 
print¢(*Enter walue for alrtenmp An"); 
ecant(rke™  belrtemp yy 
printt(PEnter value for windspd An"); 
scant (*%", dwindspd 5 
printt(*Enter value for green AW); 
scant ( "Rt", bureen y 
TREE: 


Pratt CPN nn n\n\nlnininnynlnn\nin\aiainia®); 


MODE 1; 

PrIARECMNODE WI\n"); 

{t{relhum <= 6, PS00000000E401) gato HODE2; 
else goto NODEZS; 


WODE?: 
primer (“MODE #2\n"); 
{f(pr30 «= 4, B200000000E+02) goto MODES; 
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COTA S 
else goto NOOE19; cCT.C 2 OF 7 


WODE3: - 
peiarfC“NQ0E #3\n"); 

{f(range <= 3,9000000000E+02) goto NODES; 
else goto NOOE7; 


WODES: 
printt(*NODE #6\n"); 

if(soler <= 3.0900000000E+02) goto NODES; 
else gcto TMODES; 


NODES; 
princfC@WOCE #S\n"): 

{f(veget <= 8. 1500000000E+01) goto NODES; 
else goto TNODES; ; 


WODES: 

Print#("NODE #6\n"); 

if(rethim <x 5.9500000000E+01) goto TNODE1; 
else goto TNOOE2; 


WOOE?: 

printf("NOOE #7\n"); 

if(windspd <= 2. 3500000000€+00) goto NODES; 
else goto NODEIL: 


NODES: 

PrintfC@NODE #8\n"); 

if(sr120 <= 1,0900000000E+02) goto NODE9; 
else goto NODETI; 


WODE9: 

PrintfC"NODE #9\N"); 

if(winddir <= 1. 0300000000E+02) goto TNODES; 
else goto NODE10; 


WODE 10: 

printf(™"NODE #10\n"); 

if(winddir <= 1,2600000000E+02) goto TNOOES; 
else goto TNOOE7; 


WODET1: 

printf#C@MODE #11\n"); 

if(aictemp <= 2,9000000000E+01) goto MODET2; 
etse goto NODES; 


MODE 12: 

peinté(PNOOE #12\n"); 

{f(winddir <= 2.9800000000E+02) goto TNODEB; 
else goto TNODE?; 2 


mOOE 13: 

printt(*MODE #13\n"); 

ff(airtemp <e 3. 1200000000E+01) goto TNODE10; 
eise goto TNODE1T; 


WODE 14: 
PrIATAC“MODE @14\n"); 
Af{sr60 <e 3. 6000000000€+02) goto NODE1S; 
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elee goto TNOOEI7; 


MODES: 

Printé(*MOOE #15\n"); 

iE(veget <= 7, 9300000000 ¢01) goto WODEI6; 
else goto TNOOEI; 


WODE 146: 

PRIAt(CMHOOE #16\A"); 

if(range <= 2,2500000000E403) goto WODEI7; 
wlee gato THOOE1S; 


WODEI?: 

PrinttCMWODE @17\n"); 

it (winddie <a 1,3900000000E+02) goto THODE12; 
else goto NODEI8; 


WODE18: 

PrINtf(“WOOE #18\n"); 

{f(elrremp <= 9.6000000000€+00) goto THODE 1S; 
elee goto THODET4; 


WODE 19: 

POIACECMNODE BI9\N*); 

if(windepd <= 2,8500000000€+00) goto MODE20; 
alse goto KODE24; 


WOOE20; 

Print (*NODE #20\n"); 

if(veget ce 7,.3100000000E+01) goto THOOEI8: 
else goto NOOEZ1; 


MODE 21: 

BRIACECHNODE #21\n"); 

tf (veget <= 6, 3SQ0000000E+01) goto NODEZ2; 
else goto THODE?Z2; 


MODER2: 

Printt(*WOOE #Z2\n*); 

{f(alrtemp <= 3.0900000000€+01) goto NODE23; 
else goto TNODE21; 


MODE23: 

PrINttCMHODE #23\n"); 

if(soler «# 5.6300000000€+02) goto THODEI9; 
lye goto TMODE2O; 


WODE 24; 

printhCPMQDE #246\n"); 

{f(veget «#9. S1000000008+01) goto NODEZS; 
aise goto NODE27; 


MODE 25: 

PPIATEC@MODE M25\n"); 

if(erd0 <e 7. 0500000000E +02) goto NODEZ6; 
elee goto THODE25; 


WODE26: 
printf("WODE #26\n"); 
\4(solar «= 7, Q000000000€ +02) gota TNODEZS; 


cc TAU d 
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cota o4 
else goto TNOOE24; 
ccT.c 4 OF 7 


WODE27: 

printéC*MQDE #27\n"); 

if(winddir <= 2.5200000000E+02) goto TWODE26; 
else goto TNOOE27; 


WODE28: 

printf (“Node #2B\n"); 

if(veget <= 9,3100000000E+01) goto WODE29; 
else goto NOOESS; 


NODE29: 

PrintfC"NODE #29\n"); 

if(sr120 <= 1,0000000000E+00) goto WODE30; 
else goto NOOESS; 


WOOE3O: 

printfC™NOOE #30\n"); 

if(windspd <= 5.Q000000000E+00) goto NODES1; 
else goto NODES; 


*  WODE31: 
printfCMMOoE #31\n"); 
if(salar <= 3,.S5OQ000000E+00) goto MODES2; 
else goto NOOES3; 


HODE32: 

PrintfC"NOOE #32\n"); 

if(airtemo <2 6, 7000000000E*00) goto TNODE28; 
else goto TNODE29; 


WOOET3 : 

printf("NODE #33\n"); 

if(veget <= 7.5600000000E+01) goto TNOOESO; 
else goto TNODE31; 


WODE34 : 

printéC*NODE W34\n"); 

if(range <= 2.5S00000000E+02) goto NOOES5; 
else goto NOOES6; 


WOOE3S : 

printf(“NODE #35\n"); 

if(windspd <= 2.0500000000E+00) goto TNODE32; 
else goto TNODE33; 


WODE3S: 

printf("WODE #36\n"); 

{f(winddir <e 5.SSQQ000000E+01) goto TNODES; 
else goto NODES7; 


WODE37: 

printf(SNQDE #37\n"); 

if(veget <2 7.3100000000E*01) goto TNODE3S; 
else goto TNODES4; 


WOGE 3S: 
princfC"MODE W#38\n"); 
if(range <= &,1000000000E+02) goto MODES9; 


E18 
Appendix E Listing of PC-Based Program 


co TALL 


elee goto NODES); 


ccT.c & OF 7 
WODES9: 
PrIAtACMMODE #39\N"); 
{f(er60 <= 1, 8900000000E +02) goto TMOOL 37; 
else goto TNQDESA; 


WODE4O; 

printf(MNODE #40\n"); 

if(winddie <2 7, S00000000NF +01) goto THODES9; 
alee goto NODES; 


WODES1: 

printfC@MQ0E #41\n"); 

1#(sr60 <= 1,000000000E+00) goto THODE4O; 
lee goto MODES2; 


WODE se: 

Printt(AHODE #42\n"); 

4¢(aglar <= 1,0000000000E +00) gota NOOE4$; 
else goto THODES4; 


WODEAS; 

PriNt(CMMOOE W43\n"); 

if(relhum <s 8, 7000000008401) goto NODES; 
else goto TNODESS; 


WODE SS : 

PrIMThCMMODE M44\n"); 

i#(veget <@ 8. 1900000000601) gato TNODES1; 
elaa goto THODES2; 


WODESS: 

printt(*Wade #45\n"); 

if(range «2 2,$500000000C +02) gota NODES; 
alee goto NODES7; 


NODE SS: 

printt(*NODE #4d\n"); 

it(witedspd <= 9.0000009000E +00) goto TNODESS; 
@lee goto TNODESS; 


WODES?: 

PRIACEC*MODE B47 \n4); 

if(uinddir s¢ 3, S000000000E+02) goto WODESB; 
else goto TNODES!; 


WODESS: 

printt(*WODE #68\n"); 

(f(veget <= 9. 7DO00D0000E+01) goto NCDEGI; 
elee goto THODESO; 


MODE 49: 

print#c@MODe #49\n"); 

{f(veget <= 9.$300000000E +01) goto TNODES?; 
else goto MODESO; 


WODESO: 
printé(*WODE MS0\n*); 
i*(winddir <s 1,.$200000000F +02) goto TNODESS; 
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